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ABSTRACT 


This work provides an introduc- 
tion to nonlinear wave pheno- 
mena in the context of optical 
fibres. The nonlinear properties 
of optical fibres are introduced 
with basis on a discussion of 
the elementary processes in 
nonlinear travelling wave propa- 
gation. The basic equations 
and some characteristic featu- 
res of guided-wave nonlinear 
propagation are presented. 
implications of nonlinearity for 
optical transmission systems 
are briefly discussed. 


RESUMO 


Apresenta-se neste trabalho 
uma introdução aos efeitos 
ópticos não lineares no con- 
texto das fibras ópticas. As pro- 
priedades não lineares das 
fibras ópticas são introduzidas 
com base numa discussão dos 
processos físicos elementares 
que presidem à propagação de 
ondas no regime não linear. 
Derivam-se as equações bási- 
cas de propagação e analisam- 
-se aspectos essenciais da pro- 
pagação não linear quiada. 
Discutem-se brevemente algu- 
mas consequências da não 
linearidade para os sistemas 
ópticos de transmissão. 


NONLINEAR WAVE 
PHENOMENA 
IN OPTICAL FIBRES 


Maria Hermínia Marçal 


INTRODUCTION 


The term “Nonlinear Optics” refers to phenomena that occur due to the nonli- 
near behaviour of a medium with respect to optical fields. 


Due to its many applications in lightwave communications and optoelectro- 
nics the propagation of light through nonlinear media is an increasingly 
important topic. 

While providing an introduction to nonlinear wave phenomena in optical 
fibres, this work is simultaneously intended as an introduction to fundamental 
aspects of nonlinear behaviour in itself. 

By choosing the optical fibre as the interaction medium one major simplifica- 
tion is achieved: 

In fact, conditions for one-dimensional propagation are quite well realised in 
optical fibres, due to their guiding properties, and a simple physical picture of 
the consequences of nonlinearity is made possible. 

The work begins with an intuitive approach to nonlinearity and general dis- 
cussion of the conditions when a nonlinear description of a physical situation 
is unavoidable. 


Section 2. establishes the theoretical foundations for the following sections. 
A purely phenomenological description is used throughout and the relevant 
quantities are assumed to be experimentally measurable. 


Among a wide variety of nonlinear effects, three, namely Phase-Modulation, 
Stimulated Raman Scattering and Stimulated Brillouin Scattering, are analy- 
sed in Section 3., as an illustration and an extension of the coupled-wave 
approach to nonlinear interactions established in the previous section. These 
effects are particularly important in optical transmission systems since they 
can be observed in optical fibres at power levels of a few mW. 

Section 4. discusses the propagation of light pulses in optical fibres empha- 
sizing the aspects essential to the understanding of the soliton regime and 
the properties and generation of the fundamental soliton. 


The last section briefly summarizes the effects of nonlinearity in present and 
future optical transmission systems. 


1. NONLINEARITY IN TRAVELLING WAVE PROPAGATION 


In the absence of dispersion and loss any disturbance travelling in the posi- 
tive Z-direction has an amplitude y which is a function of Z-vt, where vis the 
signal speed. 

In many “slightly” nonlinear systems v depends weakly on 1, i.e.. 
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E = v=vy(l+by) by|<<l (1.1) 
cr 


The differential equation whose characteristic lines are given by (1.1) has the 


form 


Dvd bye -0 (1.2) 


= 


É JE 


which Main, [1], designates by “travelling wave equation”. 


Equation (1.2) is a nonlinear equation since the product (W xdy/d:) is invol- 
ved. For b=0 or when by is small enough to be neglected, |by|<<1, it redu- 
ces to the linear travelling wave equation 


E o Mo. 


do (1.3) 


Not all nonlinear systems are described by such simple equations but, follo- 
wing Main, we shall use (1.1) and (1.2) to show the role of nonlinearity in 
wave propagation. 

When different points on the profile of a disturbance propagate with different 
velocities, the waveform does not maintain its shape as it propagates, 
Fig. 1., becoming progressively more distorted. 


Fig. 1 Evolution of a wave with a velocity-dependent amplitude 


A wave for which |by| is initially <<1l may propagate for some distance L, 
before any distortion becomes noticeable. For L<L, the linear representation 
is therefore suitable for treating the problem. 


Since y in a waveform is always different from zero, the distance L, for 
which the linear regime no longer applies depends on the nonlinearity para- 
meter b, which “describes” the physics of the medium, as well as on the dis- 
turbance “amplitude” y. Therefore, a nondissipative medium, where y Is not 
limited by losses, is potentially a nonlinear medium. 


* “Ho x 


Another common source of distortion is dispersion, an effect which is due to 
a dependence on frequency of the phase velocity, v, (ww). One particular 
aspect, however, distinguishes nonlinearity from dispersion: sinusoidal travel- 
ling waves are immune to dispersion but not to nonlinearity. In other words, 
sinusoidal waves maintain their form in a linear, dispersive medium, but do 
not maintain their form — they become a saw-toothed wave — in a nonlinear, 
dispersionless medium. In terms of frequency this means that in the nontli- 
near medium the original frequency spectrum is progressively enriched by 
harmonics, [2]. 

In a dispersionless medium all harmonics travel with the same phase velo- 
city, remaining in phase over long distances. This “phase matching” condi- 
tion is essential for an effective transfer of energy between different fre- 
quency components. Consequently, a dispersionless medium of sufficiently 
large dimensions is potentially a nonlinear medium. 


The linear representation is therefore an idealisation, whose applicability 
depends not only on the product by but also, as we have just seen, on the 
dissipative and dispersive properties of the medium. 


*o do + 


It is an immediate consequence of the previous analysis that all the charac- 
teristics of modern single-mode, silica-based optical fibres, namely: 


— low loss (- 0.2 dB/km at A,=1.55 um) 

— long interaction length (- 20 km for long fibres, as explained later) 

— small core radius 'a' (- 2.5 to 5 um) and associated high field intensity 
Po/ra? for a given input power 

— excitation by narrow-linewidth, single-frequency lasers (<1 MHz for cohe- 
rent detection) 


all point out to the appearance of nonlinear phenomena at relatively small 
power levels, despite the intrinsically small nonlinearity of silica. 

The long interaction length results from the fact that the conditions for one- 
dimensional wave propagation are quite well realised in optical fibres. 
Because this aspect is, ultimately, what distinguishes the optical fibre from 
bulk silica, it is instructive to consider it with more detail, 


2. BASIC EQUATIONS OF NONLINEAR INTERACTIONS IN OPTICAL 
FIBRES 


2.1 Effects of nonlinearity on plane waves 


When considering electromagnetic wave propagation in a medium it is 
necessary to specify how the electric and magnetic inductions are related to 
the electric and magnetic fields. 

Optical fibres are nonmagnetic media (u=4,) for which Z(7,0)= u,H(F,t) . The- 
refore, the only relevant quantity to consider is the electric polarization, 
p(r.1), and the electric induction is given by 


DF )=e Er)A+ P(F,) (2.1) 
Substituting (2.1) in Maxwell's equations the wave equation for the electric 
field becomes 
“E op 


Vx(VxE — == — 
CE ear ra (2.2) 


Linear optics assumes a linear relation between ? and E, P=?,, When 
nonlinearity begins to play an important role, a nonlinear term must be inclu- 
ded, P=?P +?, , to account for nonlinear effects. We then have 


Vx(VxE)+ u, “o =-H, a (2.3) 


where D(r,)=e EM.D+P (FA). 

In the regime of weak nonlinearities (i.e., when the electronic structure of the 
medium is perturbed only slightly by the applied field) 7, may be expanded 
in powers of the electric field. However, in media with an inversion 
symmetry, such as silica, all even terms in the expansion are zero ('). Con- 


sequently the relevant contribution to Pi in silica-based fibres comes from 


(') Except under very special conditions, [3]. 
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the third-order term. Higher-order terms are much weaker and will be neglec- 
ted in the present analysis. 


+ * + 


For simplicity we initially treat the electric field as a plane wave with angular 
frequency «w, propagating in the z direction. The dependence on the trans- 
verse coordinates will be accounted for later. Assuming that the field direc- 
tion is along the x axis, £=€Te, we have 


Elztyj= +24 zjexpljmo d— jko)+ ar (2.4) 


The inclusion of the complex conjugate terms (c.c.) is necessary to account 
for many nonlinear polarization terms, as will become evident later. The 
dependence of the wave amplitude on z. E(z), accounts for wave interactions 
in the nonlinear medium. 


Neglecting terms higher than third-order we obtain 


Pprse ge 


(2.5) 


E' =1 E E (=)exp(3jo, 1-3 jkz)+e.c.]+ 4 EI = Elzjexplio = jko)+CL. ) 
which shows that ?, has a term oscillating at q, and a term oscillating at 300,. 
The third-order susceptibility 7“(mw,) is in general a 4th-rank tensor but in iso- 
tropic media the number of independent terms reduces to one. 
In what follows we will concentrate on the q, term which describes a self- 
interaction effect. The 3, term is negligible in optical fibres except in very 
strict (phase-matching) conditions, [3]. 
When (2.4) is substituted in the wave equation the 1st term becomes 

PE VIRE dE aii. e 
A — al Ea — E — K E espia, — Jkz)+ ce (2.6) 

At this point one simplifying approximation can be made: it is assumed that 
the field amplitude is a slowly-varying function of Z such that the 2nd deriva- 
tive with respect to Z is negligible: 
dE 
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dE 
dz” 


With this slowly-varying amplitude approximation and *, in the form 
Pa EE ale E DE Elzjexpljm, = jkz)+ ce (2.8) 


the wave equation reduces to the following 1st-order differential equation (2) 


dE Retrato 
o == JOE (2.9) 


where ôk is a small quantity, dk = (3x /8k)(m [c)E[. In this relation k deno- 


tes the wavevector in the absence of nonlinearity, k=VI+27, o fe=no e, 
where x, is the linear susceptibility at angular frequency q, and n, is the 
linear refractive index. 

According to (2.9) the : dependence of the electric field is 
Elz)= E(z=0Jexp(-jók:). Substitution in (2.4) shows that the effect of weak 
nonlinearity on a plane wave is a small perturbation in the wavevector which 


(*) The c.c. terms are not included since they produce in this case a separate, redun- 
dant equation, 


changes from kto k+ók. This change may be interpreted as due to a refractive 
index of the medium dependent on the wave intensity 


n=n,+n!E[ (2.10) 


the nonlinear part of which is ni =(3x"/8n JE| . Equation (2.9) then takes 
the form 


IE ar. 
E =— (0 /cmiE| E (2.11) 


da. 


2.2 The influence of a transverse variation in the field amplitude 


Optical fibres, in fact, act as guides of electromagnetic waves by confining 
the wave in the two directions transverse to the propagation direction. The 
previous assumption that the wave amplitude does not depend on the trans- 
verse coordinates is, therefore, not realistic. 


As it is well known, in the linear regime any electric field component of a gui- 
ded wave can be written in the general form 


tr | a pele ti A 
E = alto EF (x yjexplgo,— jk 2)+ec.] (2.12) 


where (i) denotes the waveguide mode. The field distribution in the trans- 
verse plane, F''(x,y) and the mode wavevector at frequency o,, k''(w,) are 
determined by substituting (2.12) into the homogeneous ( Z4 =0) wave 
equation, which then yields 


E ml 3 miii Z E) 
Cd ci [= O, =" |po=0, (2.13) 
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E 
and solving this equation with the proper boundary conditions, [4], [5]. 

The form assumed in (2.12) results from a particular feature of guided 
waves: the “rigidity” in the transverse profile of the wave. In other words, the 
transverse profile of a guided wave is maintained during propagation. 

We shall now assume that the electric field is linearly polarized (*), that the 
transverse profile has been found, and that the mode wavevector has been 
determined, for a particular mode, in the linear regime. 

Treating the nonlinearity as a small perturbation we may assume that the 
refractive index of the fibre core/cladding changes from n, to n,+ôn, with 
| ônl <<n,. 


As a result, small perturbations, 6F” and ók”, occur in the transverse profile 
and in the mode wavevector, but it can be shown, [6], that one needs not to 
know the perturbed field profile to evaluate the wavevector to Íst order in ôn. 
In the same order of approximation the electric field has the form 


lr. E EI) 
| = [E (Fx )exp(ja,t— jk )+ee. (2.14) 


Substituting (2.14) into the wave equation one obtains, in the slowly-varying 
amplitude approximation 


EO Sp mê a 
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Multiplying both sides by F and integrating in the cross-section of the fibre 
yields 
E o (2.16) 
dz É 


(*) Which is valid for weakly-guíding, straight fibres, with ideal circular cross-section. 
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where n, is a nonlinear coefficient depending on the mode transverse profile 


F'da 
ny =, FE (2.17) 
[freada 


For simplicity, in (2.16) the explicit dependence on the coordinates and the 
mode index have been suppressed. 

Equation (2.16) has the same form as the equation obtained for plane 
waves, (2.11), with n., in the place of n',. The transverse profile remains 
unaffected, which means that in the absence of losses the average intensity 
of the beam (average intensity=total power/transverse area) is maintained 
during propagation. This independence on position of the intensity of the 
beam does not occur in bulk nonlinear media. 


2.3 The influence of dispersion in slowly time-varying amplitude fields 


By superimposing solutions of the linear wave equation for different frequen- 
cies we can construct wave packets that describe the propagation of optical 
pulses in the fibre 


Elzt)= ER ca jexpl item, +82) — n(k + Q)2]+ e 
É Am (2.18) 
where k+0=k(q) is the wavevector at frequency w= mw, +14). 


Although F depends on k (mw) and therefore on the frequency, we can 
assume that k(w)=k(,) for quasi-monochromatic waves (£l<<w,) and take F 
out of the integral. We then have 


| ro. 
E=—| Fix v)E(e.t)explym t— jko)+cc. 
| i | (2.19) 
where E(z,t) is a slowly-varying envelope function, [7], 


Io - 


Expanding the wavevector to 2nd order in a Taylor series about q, yields 


Limi 
O=koj-io j=t(o j+—-DKt(o, | 
(0,) = (0) +52 (0,) ni 
with 4” and k” the first and second derivatives of k with respect to «w, calcu- 
lated at (n,. 
It is now possible to differentiate (2.20) with respect to z and t and verify that 
E(z.t) is a solution of the following 1st-order partial differential equation 
E OE PE 
— + k i—— =0 
d aa "2 (2.22) 


This equation could also have been obtained in the time domain by substitu- 
ting (2.19) and (2.20) directly in the wave equation and using the slowly-var- 
ying amplitude approximation, [8], [9]. 

2.4 The joint influence of nonlinearity, dispersion and loss 


The combined effecis of dispersion and nonlinearity (*) can be obtained by 
simple addition since we are dealing with small perturbations. 


(*) We are neglecting dispersion effects in the nonlinear polarization. 


Assuming further a power loss coefficient a, losses can be accounted for in 
the usual manner, by adding a term proportional to -(o/2) times the ampli- 
tude. Taking (2.16) and (2.22) into consideration we finally have 
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dk: dt 2 x (2.23) 


This equation describes nonlinear pulse propagation in the optical fibre, in 
the 2nd-order dispersion approximation. 


3. DOMINANT NONLINEAR EFFECTS IN OPTICAL FIBRES 


3.1 Self- and cross-phase modulation (SPM and XPM) 


These two effects are due to the nonlinear dependence of the refractive 
index on the electric field amplitude. 

HW IEf has spatial variations the refractive index will also be space depen- 
dent, increasing slightly in regions of high optical intensity when the nonli- 
near coefficient is positive. This causes self-focusing in bulk media, i.e., a 
beam with amplitude larger at the centre than at the edges decreases in dia- 
meter, as it propagates. 

In optical fibres, as shown in Section 2., the transverse beam profile is sta- 
ble. However, the temporal variation of |E| will cause the phase of the peak 
of the optical pulse to be delayed with respect to its trails during propagation 
along the fibre. This effect, self-phase modulation, can therefore be conside- 
red the temporal analogue of self-focusing, [10]. 

A light wave propagating through a single-mode fibre acquires a phase shift 


dj — =( (0. fc)ndz (3.1 ) 


Integrating over the length of the fibre and assuming a power attenuation 
constant q yields 


aó=[ dó=-(0,1c)n,L+n.E(O.0fL) (3.2) 


where L, =(I-exp(-oL))/a is an effective length which replaces the actual 
length L for fibres with appreciable loss. 

The phase shift is sometimes expressed in terms of the average intensity 
I=P/A, . Defining the mode effective area A, as As =||F'da=<F' > 


and using an average refractive index n = [[nF'dajA, yields, for the tem- 
porally-varying part of the phase (º) 

Adr) = — Se N, > Exa 
“a (3.3) 
where N, is the nonlinear coefficient expressed in m?/W. N, and the previous 
coefficients n', and n, are related by (see eq. (2.17)) 


N,=2Zni/n=(2Z,/0)n, (3.4) 


(*) We are using the notation |iF"da=<r" >. The average wave intensity is 
I=P/A,=|E/n/22 [17]. 
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This phase shift causes an instantaneous frequency shift at any point in the 
pulse given by 
dAQ(!) 

dt [30] 


AD=0-W = 


which, as might be expected, is proportional to the signal power, the nonli- 
near coefficient, the effective length and the inverse of the mode effective 
area. This area is approximately the core cross section for the fundamental 
mode at a normalized frequency V=2. 


+ + * 
The theory of the previous section can be extended to the case of two opti- 


cal pulses propagating in a single-mode fibre. We shall consider only the 
case of two co-propagating waves of different wavelengths: 


E=+F(x,y | E, expljm,— jk2)+ E, expljo,t — jkazj+es .) (3.6) 


The same transverse profile has been assumed for both waves. 
The nonlinear contribution to the refractive index comes from the terms of 
the nonlinear polarization oscillating at frequencies q, (j=1,2). It is given by 
PsteE, Fa) E, “42 E 1 exp(jwm— jk z)+c.e.] 
= EX E | +2E, ) [E,exp € (87) 
The index change An, induced by the third-order nonlinearity is obtained as in 
Section 2.2 
A = n|E,| + e. 
da | (3.8) 


where n,, given by (2.17) has been assumed to be the same for both waves. 


The time-dependent part of the phase-shift as wave / propagates inside the 
fibre is now 


EM, . 
Aq (t) = = nllE +REs, 
R 


Je (3.9) 


In addition to the SPM effect, the first term in 4ó,(1), we have now a second 
term which results from phase modulation of one wave by another, co-propa- 
gating, wave. This crossed effect, cross-phase modulation, is twice as effec- 
tive as self-phase modulation for the same wave intensity. 

The propagation equations for the two optical fields 


E, .,0E o KoTE o 
= SA E | a 
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are obtained as in 2. and under the same assumptions. 


3.2 Stimulated scattering processes 


Up to now we have restrained to the analysis of nonlinear polarization terms 
which vibrate with the same frequency as the incident field. 


A second class of effects is described in the following, in which power laun- 
ched into the fibre at wavelength A, is transferred to a set of longer wave- 
lengths determined by certain characteristic vibrations of the fibre materials. 
In other words, the fibre acts as an active medium, imposing its own charac- 
teristic frequencies «, and giving origin to scattered radiation at the shifted 
frequency w,-m,. An up-shifted wave, w,+0, is also possible but that is not 
the dominant process in optical fibres and will not be considered here. 


This second class of effects includes Stimulated Raman Scattering (SRS) 
and Stimulated Brillouin Scattering (SBS), which can be classically viewed 
as a coupled three-wave interaction among an incident wave (Pump wave), 
a signal wave (Stokes wave) and a molecular vibration (Raman effect) or 
acoustical wave (Brillouin effect). As a result, part of the energy of the pump 
is progressively transferred to the signal wave. 


If no signal wave is deliberately injected into the fibre, a weak signal due to 
spontaneous emission is always present and can be significantly amplified, 
at the expenses of the pump wave. 

In order to analyse a stimulated scattering process we start with the differen- 
tial equations that describe the coupled effects between the pump and signal 
waves in the optical fibre. With 


| | | F | 
E = Fls YE, exp(jm,t— jez)+c.) E, =>F (x, YE, exp(jo,t- jk,2)+ee. (3.12) 
denoting the signal and pump waves respectively we have, in steady-state 
(d/dt=0) and for a forward scattering process 
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e en (3.13) 
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F—t=-j EH O ap EF 
dz E 8n 


where the refractive index and the transverse profile were assumed to be the 
same for both frequencies. y! denotes the susceptibility associated with the 
process. 

In terms of wave intensities /=nlE[/(2Z,) the coupled equations take the 
following, more convenient, form 


da gl 1-0, 
EA (3.14) 
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where g is a gain coefficient which is related to the imaginary part of y'*, and is 
an experimentally measurable quantity The conservation relation 
dil,+!,)ldz=0, valid in the absence of loss, was used in the derivation of 
(3.14). Subsequentily two terms «,!, and «,!, have been added to account for 
linear losses. 


The initial growth of the Stokes wave can be calculated neglecting pump 
depletion. We then have [,(:)=1,exp(-0,2), 1,=1,(2=0) and substitution in 
(3.14) vields 

IKL)= 1 (0)explgl,L —0,L) (3.15) 
where Li =(I-exp(-a,L))/a, as before. 
Equation (3.15) shows that the Stokes wave has an initial exponential 
growth, which depends on the gain coefficient, g, but also on the interaction 
length, L,;, and on the intensity of the pump, which is approximately the 
power P, divided by the effective mode area, (º) as defined in Section 3.1. In 
low-loss, single-mode fibres, the power and gain coefficient necessary for 
significant Stokes amplification can therefore be considerably reduced, com- 
pared to what happens in bulk media. 


(º) In fact we have P/A,=I=min. 
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Equation (3.15) assumes that 1,20 in z=0. However in the absence of any 
input signal the Stokes wave grows from spontaneous scattering. If the 
effective length and pump power are large enough, the Stokes wave can 
attain a magnitude comparable to the pump wave, initiating the stimulated 
process. Because the signal amplification obeys the exponential law, there is 
only a small difference in pump power between negligible conversion and 
near total conversion of pump to signal and it is possible to define a thres- 
hold or critical power, [11], defined as the input power which leads to equal 
signal and pump powers at the fibre output end. 


A similar analysis can be carried for a backward scattering process, in which 
case the following equation 


ai = =81 , +0 dl. 


dz (3.16) 
describes the evolution of the Stokes wave propagating in the negative Z 
direction. 
The critical power for stimulated forward Raman scattering has been estima- 
ted theoretically and is given by 


E = L6DA,, [EL (3.17) 


which shows a good agreement with experimental results [12]. The factor D, 
t<b<2, takes into account the relative polarization between the signal and 
pump waves: b=1 if the two waves have the same polarization, b=2 if the 
relative polarization is not maintained. Using typical values for silica fibres at 
Ao=1.55 um, q, = 0.2 dB/km, long fibres for which Lay= 1/a, = 20 km, an 
effective area A,,= Core Area=50 um” and go = 7x10"* m/W yields a critical 
power for forward Raman scattering P.=600 mW, for b=1. 
The critical power for stimulated backward Raman scattering is about 25% 
higher, which is why backward SRS is generally not observed in optical 
fibres. 

* Es * 
A strong optical field causes changes in material density through the process 
of electrostriction. A material density wave then results that propagates at 
the velocity of sound in the medium, v,. Since the refractive index changes 
with material density, the acoustic wave acts as a moving refractive index 
grating. This structure diffracts a portion of the pump wave. The scaitered 
light (the Stokes wave) is Dôppler-shifted in frequency with a frequency shift 
proportional to v,. The process is known as Brillouin scattering. 


The frequencies of the three waves are related by 
M=0-0, (3.18) 
and efficient interaction occurs under phase-matched conditions 


da =ko=k, (3.19) 


Because acoustic frequencies are much lower than optical frequencies, 
Q<<m,, O, this relation reduces to 
da = K kép — 85) (3.20) 


where e,.e, are unit vectors in the direction of the pump and Stokes waves 


respectively. Denoting by 6 the angle between the two, c,:c, =cos6 , yields 


q, = 2kpsinb/2 (3.21) 


In optical fibres the only possibilities are 9=0 and 6=7. Equation (3.21) shows 
that Brillouin amplification occurs only in the backward direction. (") 


1, (0) = É, (L) 


Op 
é E j 
I, (0) 1; (L) 
1, (0) I, (L) 
O, 


Elas 


Il, (0) I, (L) 


Acoustic Wave 


Fig. 2 Geometry of: Stimulated Forward Raman Scattering (above) and Stimulated Back- 
ward Brillouin Scattering (below) 


The evolution of the Stokes and pump waves under steady-state conditions 
are described by (3.16) and (3.14b) with 9g=g, denoting the Brillouin gain 
coefficient. Due to the small value of the Brillouin shift two simplifications can 
be made: ,/w,=1 and a,=0,=0. The initial growth of the Stokes waves is 
obtained from the resulting coupled-intensity equations in the assumption of 
negligible pump depletion. It satisfies the relation 


140)= 1 (L)explgal, Ly — OL) (3.22) 


A critical power for SBS can be estimated in a similar manner to the SRS 
case. It is found to be 


P=21hA,, | Egb (3.23) 


where the intervening quantities are defined as in (3.17). Using the peak Bril- 
louin gain for silica fibres, gs = 5x10'' m/W, and the same parameters as 
before we have P.,= 1 mW when D=1. 


The consequences of SRS and SBS for optical fibre transmission links are 
briefly discussed in Section 5. 


4. PULSE PROPAGATION IN OPTICAL FIBRES 


In describing pulse propagation in the optical fibre we shall use the equation 
derived in Section 2. for the slowly-varying envelope. This equation assumes 
that the electric field has the form of (2.19) which we now rewrite taking 
explicit account of the linear losses in the medium 


= [Ft x, v)E(z,t)exp(-oz /2)expíjm t— jkz)+ cu] (4.1) 
In the moving reference frame that follows the pulse as it propagates with 
the group velocity v,=1/k', equation (2.23) takes the form 


dE Gs KT E o ) 
E JP fa li AR fe DOE (4.2) 
= ata a do A 


(7) A very small amount of forward Brillouin scattering has been observed in optical 
fibres,[13], which is due to a relaxation in wavevector selection rules. 
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The new variable T=t-z/v, denotes the time parameter in the moving frame. 


The wave envelope can be normalized and related to the peak input power 
in the pulse, P,. We then have, using the mode effective area, Asp 


+] = — [75 74 - = itF= Eq 2] . 
E(z,7)=22,P, KnA,) expl-oz /2)u(z,T) (4.3) 
where u(z,T) is a normalized field function that includes all the shape and 
phase information in E(z,T) and n is the average refractive index. Substitu- 
tion in (4.2) yields 


7 sa 

ad IyP expl-oe jul u (4.4) 
where y is a nonlinearity coefficient, y=(m,/cKN,/A,), expressed in the 
mode effective area and the nonlinear coefficient No. 
Equation (4.4) can be reduced to dimensionless form by introducing a nontli- 
near length parameter Ly=1/(yPo) and performing the following change of 
variables: 
) t=TTo(tz/v,)/ To where To is the initial pulse width; 
li) É = z/L,, where Lo=T,*/ |k”| is the dispersion length. 


We thus obtain 


L ; a 
— + — = —— — — exp(-oz Ju) u = 0 
o dE E Ike” dt Es t pt ) | é (4 5) 


In the special case w=0 equation (4.5) is referred to in the literature as the 
nonlinear Schrôdinger equation, [18]. 


The ratio R'=L,/L,, provides a measure of the relative strengths of the nonli- 
near and the dispersion effects. The dispersion-dominant regime is characteri- 
zed by R'<<1, and using typical values for silica fibres, y =2.7 W'km' for an 
effective mode area A,,=50 umí and 4,=1550 nm this condition is easily satis- 
fied for pulses with peak power 0.1 W and initial width T, equal to 1 ps. 

The nonlinearity-dominant regime is applicable when Rº>>1, which is easily 
satisfied for pulses of initial width 10 ps with peak power of 1 W. 

The evolution of a pulse in the nonlinearity-dominant regime is described by 


E - exp(=cz)Jul u 

“e | (4.6) 
whose solution 

uz. =u(0,T)exp(jô,) (4.7) 


shows that the intensity and the pulse shape in the time domain remain 
unchanged during propagation. 


The intensity-dependent phase shift, due to SPM 


Ou = —u(O, r UE ce / Ly ) (4.8) 


is linear on the distance Z for negligible loss, az<<1 


Zm=(I-expl-0z))/a =: (4.9) 
For non-negligible loss z,y plays the role of an effective distance, smaller 
than Z. The time derivative of dy yields the instantaneous frequency shift 
from the central value q,: 


Pow AuO,T) Zg 
dl 


º oT ly (4.10) 


This equation shows that the instantaneous frequency changes across the 
pulse, i.e., the transmitted pulse is phase-modulated. A red-shift occurs at 
the leading edge of the pulse (T<0) and a blue-shift occurs at the trailing 
edge of the pulse (T>0). 


Fig. 3 Evolution of an optical Gaussian pulse in the regime of pure nonlinearity showing, 
enhanced, the instantaneous frequency shift in the output pulse (below). 


Also, the shift increases linearly with Z,, causing a spectral broadening that 
increases with distance. This is exactly the qualitative behaviour expected 
from a wave in a pure, nonlinear medium, as stated in Section 1. 


do od 


The combined action of dispersion (K"=0) and SPM is described by (4.5). Due to 
SPM, the instantaneous frequency varies across the pulse; on the other hand, 
due to dispersion, different frequency components travel at different velocities. 


In the regime of normal dispersion (k">0) lower frequencies travel faster. 
The leading edge of the pulse (red-shifted) travels faster than the trailing 
edge (blue-shifted), causing a temporal broadening of the pulse. 

Inversely, in the regime of anomalous dispersion (k"<0), lower frequencies 
travel slower, and the temporal width of the pulse becomes smaller. This 
narrowing effect is opposite to the pulse broadening effect which would 
result from dispersion alone (º). 

It may happen that the pulse amplitude, shape and central frequency are such 
that the linear broadening is exactly cancelled by the narrowing effect due to 
nonlinearity. In such case the pulse would propagate without change in shape. 
Such solution exists for equation (4.5) in the absence of loss (a=0) and is 
called the fundamental soliton. 


4.1 The soliton regime. Fundamental solitons. 


For a=0, k"<0, equation (4.5) takes the form 


dC i 207º (4.11) 
The fundamental soliton solution of (4.11) is obtained under the condition 
R=1. The corresponding profile, u, must have the form 


u=V(r)exp(-jó(6)) (4.12) 


(*) Depending on the initial pulse chirp, an initial harrowing effect may occur, [10], but 
for large distances a pulse broadening is observed for both k">0 and k'<0. 
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since it describes a pulse whose shape and spectrum are invariant over time 
and distance. It is precisely this stable form of propagation of the fundamen- 
tal soliton that makes it attractive for long distance transmission systems. 


Fig. 4 A secant-hyperbolic soliton in the z, t plane 


The amplitude V and phase à of the fundamental soliton can be obtained by 
simple integration and using adequate boundary and normalizing 
conditions (*), yielding 


u=sech(t)exp(-j0/2) (4.13) 


The peak power required for propagation of a fundamental soliton is obtai- 
ned from the condition R=1 and is given by 


E = add 


7) (4.14) 


showing that the product of the pulse peak amplitude, E, = VP, and the pulse 
width T, is a constant: a shorter pulse requires a greater amplitude to propa- 
gate as a fundamental soliton. In other words the soliton area, defined as 


Area= |  Eseen(r/T, Jar =ETA AA) 


is a constant. 

Typical values for silica at 1,=1550 nm, k'=-20 ps“/km, 7=2.7 W'km”, yield 
P;=74 mW for a pulse with T,=10 ps. 

The fundamental soliton is a special solution of (4.11). With generality the 
evolution of a pulse, u(Z,T), must be determined by solving (4.11) subject to 
the boundary condition that at Z=0 the solution must match the input pulse. 
The evolution of pulses of initial secant-hyperbolic shape but variable peak 
power has been studied in detail both theoretically and numerically [14],[15]. 
For peak power equal to NºP,, with N an integer and P, given by (4.14), the 
pulses still exhibit pure soliton behaviour but their structure changes periodi- 
cally with z (these are called N-order solitons). 

For N=1+e, le|<0.5 the pulses were shown to evolve asymptotically to a fun- 
damental soliton 


u=% sech(xT)exp(-jx'€/2) (4.16) 


(º) Vand av/dr vanish when C5=s;vV=1 and dV/dr=0 for 1=0. 


where x=I+2e is a scale factor. Part of the pulse energy is dispersed away 
in the process. 

A similar behaviour occurs for N=1 and initial shape different from the 
secant-hyperbolic profile: the pulse adjusts its shape to evolve into a soliton, 
[3]. 

The soliton is therefore a stable formation with respect to small perturbations 
either in amplitude or in shape, and can be easily excited if the initial peak 
power exceeds a threshold value. 


The effect of loss is to reduce the amplitude as the pulse propagates. If the 
fibre loss is too large the amplitude of the pulse decreases as it propagates 
and the pulse may be dissipated before it forms a soliton. The minimum 
input power required to excite a soliton is therefore aggravated by fibre loss. 


5. EFFECTS OF NONLINEARITY IN OPTICAL TRANSMISSION: CON- 
CLUDING REMARKS 


In the preceding sections we have seen that the optical fibre is an inherently 
nonlinear medium. Each of the nonlinearities described affects optical trans- 
mission in different ways. 


The influence of nonlinear effects can be both negative and positive: Severe 
degradation in a transmission system capacity can result ftom nonlinear los- 
ses, pulse broadening and crosstalk; On the other hand, when nonlinear pro- 
cesses result in signal amplification and dispersion compensation they can 
significantly increase the system performance. 


The analysis presented here is very brief, and based only on the critical 
powers ('º) estimated in Section 3. 


We can conclude that SRS is not likely to occur in single-channel optical 
communication systems, since the power launched into the fibre is typically 
around 1 mW. 

However the situation is quite different in wavelength-division multiplexed 
(WDM) systems, because numerous wavelengths are injected into the fibre 
and the longer wavelentghs will be amplified by the shorter wavelengths. 
Imposing an amplification of one channel by another less than 1% (20 dB 
crosstalk) yields g,L,P./bA, =Ol, [11], which is a condition significantly 
more restrictive than (3.17). 


For the same parameters as before and assuming maximum Raman gain 
(the worst case), a critical power of 600 mWx0.01/16=0.4 mW is'obtained. 
Although channel separation equal to the precise value that corresponds to 
the maximum Raman gain can be avoided, the large bandwidth of the 
Raman gain (=10 THz) makes this effect a major concern in WDM systems. 
The most likely nonlinear loss mechanism in single-channel optical fibres is 
SBS, for which a critical power as low as 1 mW was obtained in Section 3. 
However, maximum Brillouin gain occurs for pump lasers with linewidths less 
than the Brillouin linewidth (Af,<100 MHz in silica fibres). For lasers with line- 
widths Af,>>Af,, the Brillouin gain decreases as the ratio Afg/Af,, resulting in 
an increase in critical power. Also, due to the small frequency shift, SBS is 
not as dangerous for WDM systems as SRS, and the critical power for SBS 
is almost independent of the number of channels, [12]. 


('º) Higher-order nonlinear and dispersion effects and many other factors, including 
modulation schemes or channel separation, were not discussed. 
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The intensity-dependent refractive index origins pulse broadening and can 
be a source of noise in phase-modulated systems. 


Although no threshold power exists for SPM in the same sense as in the 
other effects, a critical power can be estimated by imposing that in a disper- 
sion-limited link, nonlinear broadening does not exceed linear broadening, 
Stolen [11]. Alternatively, for coherent detection systems where phase stabi- 
lity is important, the critical power can be estimated by imposing qy;<<1. For 
long fibres L.1/a and from (3.3) one obtains P,=o/y. Using as parameters 
q=0.2 dB/km, y = 2.7 Wº km yields P,, =17 mW, [16]. The same order of 
magnitude was obtained by Stolen, for Gaussian pulses. 

Phase noise becomes worse with WDM systems because XPM is twice 
more effective than SPM for the same amount of power. 


do * 


Among the potentially useful phenomena that arise from nonlinearity we can 
include Raman scattering and soliton propagation. 

The Raman effect is frequently used for optical amplification, due to its broad 
gain-bandwidth. 

The propagation characteristics of the fundamental soliton make it potentially 
interesting for long distance, high-bit rate optical transmission links ("). 
However, since soliton propagation is critically dependent on the intensity of 
the pulse, a soliton transmission system requires that the pulse be reshaped. 
The possibility of using the Raman effect to compensate for fibre loss has 
been extensively studied, [19], and is referred here as an example of double 
application of nonlinear propagation (soliton propagation and the Raman 
effect). Some difficulties in the practical realizability of such systems has 
recently made the attention turn to other soliton regeneration methods, [20]. 
Soliton propagation offers one solution in the future development of large- 
capacity, long-distance optical transmission. 

Other solution is WDM coherent communication systems. 

Both are the subject of great current interest, [21], and they are closely rela- 
ted since the nonlinear effects can influence strongly the coherent systems 
characteristics and vice-versa. 

In either case, nonlinear optical phenomena will have an increasingly impor- 
tant role in optical fibre communications and associated technology. 
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RESUMO 


Os resultados obtidos através 
de uma técnica óptica de visua- 
lização por plano de luz laser 
permitem analisar o efeito da 
variação da pressão de atomi- 
zação e da existência de um 
escoamento de ar coaxial na 
qualidade de atomização e no 
ângulo de sprays provenientes 
de diferentes tipos de injectores 
de pressão, nomeadamente, de 
cone oco e cone sólido. São 
estudados sprays provenientes 
de injectores de pressão na 
presença e ausência de um 
escoamento de ar coaxial, São 
elaborados vários critérios para 
determinação do ângulo de um 
spray sendo também analisada 
a incerteza associada aos 
resultados obtidos. 


ABSTRACT 


A conventional laser light sheet 
technique was used to study 
overall external characteristics 
of different sprays. Visualization 
was carried out for two distinct 
pressure nozzles (full and hol- 
low cone), in the presence and 
absence of coaxial air flow. As 
the main objective is to deter- 
mine spray angles, several cri- 
teria were used, and some 
were established based on the 
observed behaviour of the spe- 
cific flow configuration under 
study. An extensive analysis 
was performed with the objec- 
tive of establishing the most 
adequate criteria that seemed 
to be suitable for these and 
other flow configuration. 


Palavras-chave: injectores de pressão; ângulos; visualização 
Keywords: pressure nozzles; spray angles, visualisation 
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1. INTRODUÇÃO 


A atomização de um líquido pode ser obtida à custa de vários tipos de injec- 
tores podendo, contudo, distinguir-se duas classes fundamentais (Faeth, 
1977; Lefebvre, 1983): os injectores de pressão e os assistidos a ar. Nos 
injectores de pressão, o liquido passa através de um orifício de dimensões 
reduzidas, sendo a pressão de injecção um factor crítico, uma vez que con- 
diciona a velocidade que é imprimida ao líquido, parâmetro determinante da 
dimensão das gotas obtidas. Existem, vários injectores de pressão, que se 
distinguem pela existência ou não de rotação do liquido à saída do orifício. 
Os injectores de pressão sem rotação permitem obter sprays de ângulos 
reduzidos, entre 5 e 15º (Lefebvre, 1989), em que as gotas se distribuem de 
forma aproximadamente uniforme pela totalidade do volume do spray. 
Sprays deste tipo são designados por sprays de cone sólido. 


No caso dos injectores de pressão com rotação, o líquido emerge do injector 
na forma de folha de líquido anular, que se dispersa radialmente formando 
um spray oco de forma cônica, com a maior parte das gotas concentradas na 
periferia. Os ângulos típicos destes sprays situam-se entre 30 e 180º, depen- 
dendo da aplicação (Lefebvre, 1989). Os injectores de pressão com rotação 
também permitem a obtenção de sprays de cone sólido pela utilização de um 
processo de introdução de gotas no interior do spray de cone oco. O desem- 
penho dos atomizadores de pressão tem sido estudado nos seus mais diver- 
sos aspectos por vários autores (e.g. Mao, 1991; Tokuoka et al., 1991). 


Um processo alternativo de atomização, faz uso do conceito de injectores 
assistidos a ar (ver, por ex., Lefebvre, 1989) e consiste basicamente, em 
expor um escoamento de liquido de velocidade reduzida, a um escoamento 
de ar de velocidade elevada, com o objectivo de melhorar a qualidade de 
atomização. O atomizador em estudo neste trabalho pode enquadrar-se 
neste contexto, uma vez que consiste num injector de pressão com rotação 
sujeito à acção de um escoamento coaxial de ar. 


Na maior parte das aplicações a função do atomizador não é apenas a de 
promover a desintegração do líquido em gotas de reduzidas dimensões. 
sendo, também, importante que as “descarregue” no meio gasoso envol- 
vente, de forma axissimétrica. Esta característica conjuntamente com o 
angulo do spray é de grande importância nos processos que envolvem com- 
bustão, dada a sua influência no processo de ignição, nos limites de flama- 
bilidade, e na emissão de hidrocarbonetos não queimados e fuligem. 

O ângulo de um spray é um parâmetro que, por inspecção visual, tem uma 
definição intuitiva imediata, sendo determinado pela zona de fronteira entre 
o escoamento de líquido proveniente do atomizador e o meio ambiente. 
Contudo, esta definição é demasiado vaga, e pouco rigorosa, para ser utili- 
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zada como critério na determinação de ângulos de sprays a partir de resulta- 
dos de visualização. A maior dificuldade existente na medição do ângulo de 
um spray reside no facto das fronteiras do escoamento serem curvas, o que 
se deve a efeitos de interacção entre o spray e o ar ambiente. Consequente- 
mente, o ângulo do cone é dado pelo ângulo formado por duas rectas que se 
cruzam na saída do injector e cortam as fronteiras do spray a uma dada dis- 
tância a jusante da saída do atomizador (Lefebvre, 1989). A definição desta 
distância é, por vezes, função do contexto específico em que se pretende 
determinar o ângulo de um spray. É natural que se procure uma definição 
que dê maior relevância às fronteiras do escoamento junto à saída do injec- 
tor (uma vez que é aí que todos os fenómenos que determinam a qualidade 
de atomização têm lugar). No caso de aplicações em que é necessário quan- 
tificar o grau de penetração de um escoamento, assim como a sua taxa de 
abertura, o critério de determinação do ângulo podera basear-se nos contor- 
nos do escoamento em regiões mais a jusante da saida. 

No trabalho aqui apresentado é dispensada uma atenção especial à elabo- 
ração de critérios de determinação de ângulos de sprays com base em dois 
tipos de sprays: os sprays de cone sólido e os de cone oco. Este artigo 
encontra-se dividido em 4 partes, incluindo esta introdução. A secção 
seguinte refere a instalação experimental e o atomizador usado, conjunta- 
mente com a técnica experimental utilizada na caracterização dos escoa- 
mentos a jusante do atomizador. A secção 3 discute a definição dos critérios 
de determinação de ângulos e sua aplicação aos escoamentos estudados. 
Nesta secção está também incluída uma análise da influência do critério 
usado nos resultados obtidos, que possibilita a elaboração de algumas con- 
clusões sobre a determinação de ângulos de um spray. A última parte 
resume as principais conclusões do presente estudo e sugere algumas 
áreas de interesse para trabalho futuro. 


2. INSTALAÇÃO EXPERIMENTAL E INSTRUMENTAÇÃO 


2.1. Instalação Experimental 


A instalação experimental utilizada foi concebida de modo a possibilitar o 
estudo de diferentes tipos de atomizadores comerciais, assim como de 
modelos laboratoriais (Carvalho, 1995). O atomizador estudado encontra-se 
representado na figura 1, sendo constituído por dois cilindros concêntricos. 
A parte central, na qual se incorpora o injector de pressão com rotação, tem 
um diâmetro exterior, D,, de 12.0 mm e a parte exterior correspondente ao 
escoamento de ar tem um diâmetro externo, D,,, de 26.7 mm. 


Liquido 


Ar coaxial 


D,=26.7 mm 


Fig. 1 Esquema do atomizador 


O sistema de alimentação de líquido utiliza água como fluido, sendo a sua 
alimentação efectuada a partir de um tanque de pressão, que permite obter 
pressões de água à entrada dos injectores entre 2 e 5 bar, com um erro 
máximo de 5%. O erro na medição do caudal de liquido à saída do injector 
é de 1%. 


Foram estudados dois injectores distintos, nomeadamente, um injector de 
cone sólido (F60), com um ângulo nominal de 60º e um caudal de líquido de 
22.90 g's, para a pressão de 5 bar, e um de cone oco (H80), com um 
ângulo nominal de 80º e um caudal de líquido de 2.77 g/s, também para 
P = 5 bar. Os resultados obtidos do procedimento de calibração destes 
injectores permitiram concluir que a pressão de 5 bar é a pressão de funcio- 
namento mais indicada, uma vez que é aquela para a qual os valores de 
caudal de líquido e de ângulo do spray mais se aproximam dos valores indi- 
cados pelo fabricante. 


2.2. Instrumentação 


A visualização dos escoamentos foi efectuada recorrendo a um plano de luz 
laser, de 2 mm de espessura, que foi obtido pelo espalhamento de um raio 
laser, com uma potência máxima de cerca de 5 W, ao passar através de 
uma lente cilíndrica. O plano de luz atravessa o escoamento segundo o eixo 
de simetria do injector. Os resultados foram obtidos recorrendo ao uso de 
uma máquina fotográfica e filmes de 400 ASA, com tempos de exposição 
entre 8 e os 125 ms e com um factor de escala de 1:2.7. O tempo de expo- 
sição foi escolhido de acordo com as características do escoamento em 
estudo. A intensidade de luz do plano de luz laser foi optimizada para cada 
caso de forma a garantir o tempo de exposição adequado. Recolheram-se 
três fotografias, com tempos de exposição distintos, para cada condição de 
funcionamento. 


3. APRESENTAÇÃO E DISCUSSÃO DOS RESULTADOS 


Nesta secção apresentam-se alguns resultados da visualização dos escoa- 
mentos estudados, de forma a permitir a elaboração e aplicação de diversos 
critérios de determinação de ângulos de diferentes sprays. Apresenta-se, 
ainda, uma análise da influência do critério usado na determinação do 
ângulo nos resultados obtidos. Os resultados aqui apresentados são refe- 
rentes aos dois injectores ensaiados: um de cone sólido e outro de cone 
oco, para pressões de líquido de 2,3, 4 e 5 bar, na presença (U,= 30 m/s) e 
na ausência (U, = 0) de ar coaxial. 


Tipo de spray - Definição do ângulo 


As figuras 2 e 3 apresentam os resultados da visualização dos sprays de 
cone sólido (injector F60) e de cone oco (injector H80), respectivamente. 
Estes resultados ilustram claramente a diferença entre o tipo de sprays em 
estudo. No caso do injector de cone sólido observa-se uma distribuição uni- 
forme de líquido ao longo da secção transversal do escoamento. No injector 
de cone oco, é na periferia do escoamento que se encontram concentrações 
elevadas de líquido, tal como se pode observar na zona imediatamente a 
jusante da saída do injector. A análise destes resultados ilustra, ainda, a 
diminuição da “qualidade de atomização” para valores reduzidos de pressão 
de líquido, sendo este comportamento independente do tipo do injector estu- 
dado. As alíneas c) das figuras 2 e 3 ilustram a influência do escoamento de 
ar coaxial nos sprays obtidos para uma pressão de 5 bar, sendo esta mais 
nitida para o injector H80 do que para o F60. 
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Fig. 2 Visualização do escoamento do injector F6O para diferentes pressões de funcionamento: 
a)P=S5bar,m =22909s:b)P=2bar,m, = 13.43095,c)P=5bar,m, = 22.90 gs e 
U, = 30 m/s 


a) b) c) 


Fig. 3 Visualização do escoamento do injector H80 para diferentes pressões de funcionamento: 
a)P=5bbar,m=2.77g's;b)P=2bar,m=1.729's;c)P=5bar,m =2.77 g'se 
U, = 30 m/s 


Os resultados para os injectores de cone sólido e cone oco demonstram uma 
clara diferença de estrutura intrinseca dos dois escoamentos obtidos, factor 
importante a ter em conta na elaboração de critérios para a determinação dos 
respectivos ângulos. Consequentemente assume-se, à partida, a necessidade 
de elaborar critérios distintos para os sprays de cone sólido e cone oco. 


O critério usado na determinação do ângulo dos sprays de cone sólido 
baseia-se na zona da saída do injector e será designado por “critério da 
saída”. Esta região é facilmente identificada pela elevada intensidade de luz 
dispersada que lhe estã associada. Contudo, a aplicação exclusiva deste 
critério pode conduzir a alguma ambiguidade na determinação do ângulo do 
spray. Esta é ultrapassada recorrendo às regiões a jusante da saída onde 
se identificam as trajectórias das gotas, não necessariamente as mais afas- 
tadas da linha central, mas sim as consideradas significativas em termos de 
escoamento. Este critério será designado por “critério das trajectórias”. Na 
figura 4 a) exemplifica-se a aplicação destes dois critérios à determinação 
do ângulo do spray obtido pelo injector F60, a uma pressão de 5 bar. 

No caso dos sprays de cone oco toma-se como origem do ângulo a região 
vizinha da saída do orifício do injector, ficando este completamente definido 
com base nos limites das zonas de maior dispersão de luz, que na zona da 
saída do injector, mais se afastam da linha central. Este critério designa-se 
também por “critério da saída”. 
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Fig. 4 Critérios de determinação de ângulos de sprays: 
a) Injector F60 (P = 5 bar; U, = 0); b) Injector H80 (P = 5 bar; U, = 30 m/s) 


A introdução de ar coaxial origina uma maior concentração de líquido perto 
da saída, devido ao facto daquele escoamento de ar condicionar fortemente 
a trajectória das gotas que passam a acumular-se numa região mais prô- 
xima da linha central (ver figura 3 c). Verificou-se que a aplicação do “critê- 
rio da saída” conduz a uma sobreestimação dos valores dos ângulos, 
essencialmente para os escoamentos com ar coaxial. Consequentemente, 
estabeleceu-se um outro critério que consiste em definir o ângulo do spray 
com origem na saída do injector na zona em que existe spray 1 D.y a 
jusante da saída. Este critério é designado por “critério 1 D.w a jusante da 
saída”. Na figura 4 b) exemplifica-se a aplicação do critério “1 D,, a jusante 
da saída” e o “critério da saída” na determinação do ângulo do spray obtido 
pelo injector H80 a uma pressão de 5 bar e com escoamento de ar coaxial. 


Ângulo do spray 


Da análise detalhada de um conjunto de imagens, como as ilustradas nas 
figuras 2 e 3 retirou-se a informação necessária à quantificação dos ângulos 
dos vários escoamentos, de acordo com os critérios atras estabelecidos. 
A figura 5 a) apresenta os resultados obtidos da aplicação do “critério da 
saída” conjugado com o “das trajectórias” ao injector F60. Os resultados 
mostram que à medida que a pressão de atomização aumenta, o ângulo do 
spray se aproxima do valor nominal de 60º. Verifica-se que com a introdu- 
ção de ar coaxial a diminuição do ângulo do spray é tanto maior quanto 
menor é a pressão de atomização sendo aquele decrescimo de 11.4%, para 
uma pressão de 2 bar, e de 2.8% para 5 bar. 

A figura 5 b) apresenta os ângulos obtidos para o injector H80, pela aplica- 
ção do “critério da saída” e de “1 D,, a jusante da saída”. Verifica-se que o 
“critério da saida” não parece ser o mais adequado ao escoamento em 
estudo, uma vez que conduz a uma evolução de q, que notoriamente não é 
real (ver também figura 4b). À inspecção visual do escoamento revela uma 
diminuição do ângulo do spray na presença de ar coaxial, independente- 
mente da pressão, o que está em contradição com os resultados obtidos 
pelo “critério da saida” para pressões de 2 e 3 bar. 


Esta situação é ultrapassada quando se recorre ao critério de “1 D.y a 
jusante da saida”. De acordo com o critério de “1 D,, a jusante da saída” a 
utilização de ar coaxial conduz a uma diminuição no ângulo de cerca de 
93% para uma pressão de 2 bar e de cerca de 11% quando a pressão é de 
3 bar. 
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Fig. 5 Representação do ângulo do spray em função da pressão do liquido: a) injector 
F60; b) injector HBO. 


Os resultados obtidos estão de acordo com os encontrados por outros inves- 
tigadores (ver, por ex., Lefebvre, 1989) na gama de pressões estudadas. 


Influência do “critério de determinação do ângulo” nos resultados obti- 
dos 


Pelo conjunto de resultados atrás apresentado torna-se clara a influência do 
critério usado nos valores dos ângulos obtidos. Esta situação revela-se parti- 
cularmente crítica no caso dos sprays de cone oco. À diferença entre os 
dois critérios usados no caso destes sprays consiste na dimensão da zona 
do escoamento que é usada para determinar o seu ângulo. No “critério da 
saída”, aquela zona é de menores dimensões do que a considerada no cri- 
tério “1 D,, a jusante da saída” (tal como se pode observar na figura 4b). 
No sentido de avaliar a influência que a distância à saída, X, tem na deter- 
minação dos ângulos, representam-se na figura 6 resultados relativos aos 
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Fig. 6 Ângulos dos sprays relativos aos injectores F60 e H80, a várias pressões e na 
ausência de ar coaxial. 


injectores F60 e H80, obtidos a partir da taxa de abertura dos sprays em 
função da distância à saída, de acordo com a expressão: a=2-arctg( RA), 
onde R é metade da largura do spray a uma distância X, da saída do injec- 
tor. Uma análise mais detalhada da influência dos diversos critérios de 
determinação dos ângulos para diferentes configurações estudadas pode 
ser encontrada em Santos (1996). 

Esta análise mostra, mais uma vez, a diferença intrínseca do tipo de spray 
em análise. Enquanto as curvas de a em função de X relativas ao injector 
de cone oco apresentam um ponto de inflexão, na região a jusante de 
X = 150 mm, o mesmo não se verifica para o injector de cone sólido. Em 
ambos os casos, a zona de maior gradiente encontra-se na vizinhança do 
injector, a montante de X = 50 mm, independentemente da pressão utili- 
zada. A figura 6 sugere, ainda, a existência de uma distância à saída a par- 
tir da qual a determinação do ângulo se torna menos dependente de X. No 
caso dos sprays de cone sólido aquela distância deve ser superior a cerca 
de 50 mm. A região dos sprays de cone oco em que o gradiente do ângulo 
é mais reduzido situa-se entre 75 e 175 mm, sendo o limite superior depen- 
dente da inflexão da curva de a em função de X. 


4. CONCLUSÕES 


Para a caracterização dos escoamentos provenientes de diferentes injecto- 
res foi usada uma técnica de visualização com plano de luz laser, cujos 
resultados permitem ilustrar as diferenças fundamentais de estrutura 
existentes entre os sprays de cone sólido e de cone oco, bem como a in- 
fluência que um escoamento de ar coaxial pode ter nas suas características 
externas. 


Partindo de uma definição clássica do ângulo de um spray, elaboraram-se 
critérios adequados à sua medição, tendo ficado clara a necessidade de 
atender às características intrinsecas dos escoamentos no estabelecimento 
daqueles critérios. Embora em ambos os tipos de sprays os critérios se 
baseiem fundamentalmente na zona da saída, no caso dos sprays de cone 
sólido é necessário ter em conta as trajectórias das gotas em regiões dis- 
tantes da saída do injector. No caso dos sprays de cone oco foram definidos 
dois critérios alternativos, conduzindo o critério “1 D,, a jusante da saida” a 
resultados mais consistentes com a inspecção visual do escoamento. 


A principal fonte de incerteza na determinação do ângulo de um spray é 
referente ao critério de medida propriamente dito, bem como ao rigor na sua 
aplicação. A utilização de vários critérios de determinação de ângulos, bem 
como o estudo da influência da distância à saída do injector e da posição da 
origem virtual do spray, nos ângulos medidos, conduziu à obtenção de uma 
dispersão de valores, a partir dos resultados da visualização de um único 
escoamento. Consequentemente, fica demonstrada a necessidade de esta- 
belecer um critério de medição claro e universal, no sentido de se poderem 
correlacionar valores de ângulos de sprays, provenientes de diferentes injec- 
tores. 
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AS PAREDES DOS 
EDIFÍCIOS ANTIGOS 


Fernando Farinha da Silva Pinho 
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1. INTRODUÇÃO 


Passaram-se já cerca de 12.000 anos desde que as primeiras casas de 
terra foram construídas em Jericó, a cidade mais antiga do mundo [8]. 


Em Portugal, terão surgido há cerca de 2.500 anos, mas deixaram de ser 
construídas como solução habitacional, característica de populações de bai- 
xos recursos económicos e de zonas de abundante terra barrenta, por volta 
da década de cinquenta, para se voltarem a fazerem experiências recentes, 
designadamente no Algarve, tendo em vista a sua implementação, ainda 
que tímida, dados os custos inerentes, à actividade turística. 


São, portanto, uma referência histórica muito importante, quer sob o ponto de 
vista social, quer sob o ponto de vista da própria evolução da capacidade 
Humana de adaptação ao meio envolvente, se tivermos presente que uma 
habitação é, antes de mais, um local onde as pessoas tentam, no mais primário 
dos estadios, estar seguras e protegidas das condições climatéricas adversas. 
Verifica-se, por outro lado que, ao longo dos tempos, as pessoas procura- 
ram criar habitações dimensionadas e construídas de acordo com os mate- 
riais mais característicos das zonas onde se pretendiam instalar, como é o 
caso das já referidas construções de terra e, bem assim, das construções 
em pedra, madeira ou mistas. 


É, portanto, neste contexto que se desenvolve o presente trabalho, onde se 
tentam sintetizar as principais soluções construtivas existentes no nosso 
Pais nos últimos 200 anos, mas com maior incidência no primeiro quartel 
deste século, tendo como principal referência a Tese de Mestrado intitulada 
“SISTEMATIZAÇÃO DO ESTUDO SOBRE PAREDES DE EDIFÍCIOS ANTI- 
GOS — Análise de casos” [5], do mesmo autor do artigo, cuja publicação 
está a ser preparada, após algumas alterações pontuais entretanto efectua- 
das e onde se descrevem e caracterizam as principais soluções construtivas 
de edifícios antigos de habitação em Portugal, antes do terramoto de 11 de 
Novembro de 1755, até meados do século XX. 


2. PAREDES RESISTENTES E DE COMPARTIMENTAÇÃO 
DE EDIFÍCIOS ANTIGOS 


As paredes resistentes dos edifícios antigos distinguem-se em termos estrutu- 
rais das de compartimentação pelas suas maiores secções, que lhes proporcio- 
nam substancial melhoria de características mecânicas para resistirem aos 
esforços aciuantes e, de um modo geral, para a estabilidade global do edifício. 
De facto, tendo as paredes divisórias características geométricas inferiores 
às anteriores, desempenham fundamentalmente funções de compartimenta- 
ção, embora também contribuam para o travamento geral da estrutura do 
edifício, em função da sua organização ou distribuição espacial. Com o pas- 
sar do tempo, no entanto, estas paredes divisórias vão tendo um acréscimo 
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de solicitações, dadas as constantes e progressivas alterações das condi- 
ções de equilibrio estático, devidas ao envelhecimento e fluência dos mate- 
riais e/ou aos movimentos diferenciais das fundações, acréscimos de sobre- 
cargas e degradações várias resultantes de acções sismicas [1], a que a 
região sul-centro de Portugal é particularmente sensivel. 


3. PRINCIPAIS SOLUÇÕES DE PAREDES DE EDIFÍCIOS ANTIGOS 


Referem-se em seguida as soluções mais representativas de paredes de 
edifícios antigos, construídas em Portugal, nos últimos dois séculos. 


3.1 Paredes de cantaria e enxilharia 


Caracterizam os edificios mais nobres, como os monumentos e casas 
senhoriais, sendo a sua estabilidade garantida à custa da elevada espes- 
sura das paredes e dos perpianhos, pedras que ocupam toda a espessura 
da zona corrente da parede, e pela forma cuidada como eram feitos os 
cunhais, garantindo um adequado travamento das pedras, que podiam ainda 
ser gateadas, com materiais diversos (metal, chumbo, madeira, pedra), e 
solidarizadas entre si com argamassa de cal e areia. 


3.2 Paredes de alvenaria ordinária e mista 


Representam uma larga percentagem das soluções construtivas de edifícios 
antigos no nosso País, sendo os de maior porte (2 a 5 pisos) próprios das 
zonas urbanas, enquanto os de pequeno porte se localizam principalmente 
nas zonas rurais. 

A construção de paredes de alvenaria mista, principalmente de cantaria e 
alvenaria era vulgar. 


A cantaria aplicava-se para efeitos decorativos, ou para reforçar as zonas 
sujeitas a maiores esforços ou com maiores probabilidades de receberem 
choques. Quando aplicada para efeitos de adorno, a cantaria recebia o 
nome de pedra de forro, tendo uma espessura média de 0.08 m a 0.10 m. 
O revestimento era feito em simultâneo com a construção da parede, sendo 
cada pedra assente com argamassa e com o auxílio de peças metálicas 
(gatos de ferro), Invisíveis nos paramentos [4,6]. 


Mas, um dos mais representativos exemplos de paredes mistas (de alvenaria 
e madeira [3]) é o que corresponde à reconstrução da Baixa Pombalina em 
Lisboa, após o terramoto de 1755, impulsionada por Sebastião José de Car- 
valho e Melo, Marquês de Pombal e Ministro do Rei D. José |, caracterizadas 
por paredes constituídas por uma gaiola ou armação de madeira embebida 
no maciço de alvenaria das paredes mestras, à face do paramento interior, 
intimamente ligada pelas mãos (pequenas peças de madeira emergentes da 
gaiola, no sentido da alvenaria), mas de forma a que, no caso da ocorrência 
de um sismo o edifício se desmoronasse, mantendo íntegra a gaiola tridimen- 
sional de madeira, contraventada pelas cruzes de Sto. André, fig. 1. 


São ainda muito comuns em Lisboa as paredes mistas de alvenaria e tijolo 
[3], sendo este aplicado nas cadeias verticais, nas cadeias horizontais, nas 
guarnições, em revestimentos, etc.. 


3.3 Paredes de pedra seca 
São próprias das zonas onde a cal rareava, nomedamente em zonas rurais, 


pelo reduzido porte permitido, não superior a dois pisos (andar térreo e 
sobrado), para espessuras de 60 a 80 cm. 


Fig. 1 Modelo da Gaiola Pombalina, existente nos B. M. de Lisboa, para formação dos 
Bombeiros 


A predominância de xistos e granitos permitia a obtenção económica de 
grandes blocos, os silhares, de forma paralelipipédica, com que se cons- 
truíam as paredes em fiadas praticamente niveladas. Originavam constru- 
ções secas no Inverno e frescas no Verão, mas eram susceptíveis a alguns 
problemas estruturais, nos casos de alvenaria muito irregular ou mal 
assente, e de problemas termohigrométricos, nos casos em que contacta- 
vam permanentemente com solos permanentemente saturados, porque a 
grande capacidade de armazenamento podia tornar-se irreversível [1]. 


3.4 Paredes de alvenaria aparelhada 


A alvenaria aparelhada corresponde ao antigo opus incertum ou aparelho 
poligonal dos Romanos [6]. Era constituída por pedras irregulares aparelha- 
das numa das faces e assentes em argamassa, seleccionando-se as mais 
rijas e de melhor aspecto para os paramentos, originando assim uma super- 
fície praticamente plana. As arestas eram aperfeiçoadas, dando-lhes formas 
regulares. 


3.5 Paredes de tijolo 


O tijolo mais usado na construção dos edifícios antigos (tijolo burro ou 
furado), tinha em média as dimensões de 0.23 m x 0.11 m x 0.07 m, as jun- 
tas não excediam 0.01 m e a argamassa de aplicação para além de ter 
alguma consistência devia também preencher todos os espaços livres, 
sendo para isso chapada antes do assentamento do tijolo, que se carregava 
e esfregava levemente até a argamassa jorrar pelas juntas, [4,6]. 


No assentamento do tijolo, devia-se ter em conta se ele era assente à vista 
ou não, pois no primeiro caso a escolha devia ser muito cuidada, rejeitando- 
-se os fendidos, partidos, mal cozidos, queimados e os que tinham dimen- 
sões anormais que eram reservados para as partes da parede onde não 
ficavam à vista. Se a parede era rebocada, todos os tijolos serviam, não 
havendo sequer grande preocupação com as irregularidades das juntas que, 
aliás, contribuiam para melhorar a aderência da argamassa do reboco [4]. 
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Em tempo quente e seco, molhavam-se os tijolos para lhes tirar a poeira, 
facilitar a aderência da argamassa e evitar que absorvessem de imediato a 
sua humidade. No assentamento dos tijolos, deviam-se matar as juntas, evi- 
tando-se, como acontecia no caso das pedras, a correspondência das juntas 
verticais, existindo para tal várias formas de os dispor nas fiadas, dependen- 
tes fundamentalmente da espessura da parede. 


Em igualdade de circunstâncias, as paredes de tijolo podiam ser menos 
espessas que as de pedra, porque a horizontalidade dos leitos e a maior 
regularidade na sobreposição lhes conferiam maiores resistências; além 
disso tinham uma condutibilidade térmica inferior. Assim, uma espessura de 
0.35 m a 0.38 m numa parede de tijolo, correspondia a uma parede de alve- 
naria corrente com 0.50 m [2]. No QUADRO |, referem-se as designações 
atribuídas às paredes de tijolo, em função da sua espessura (ou seja: do 
número e disposição dos tijolos nelas incluídos). 


QUADRO | - Designação das paredes de tijolo, 
em função da sua espessura, [4,6] 


ESPESSURA DA PAREDE DESIGNAÇÃO 


Igual à espessura de um tijolo Pano de tijolo ao alto 
Igual à largura de um tijolo Parede de meia vez 
Igual ao comprimento de um tijolo Parede de uma vez 
igual à soma do comprimento com a largura E 
| de um tijolo o 
| Igual duas vezes o comprimento de um tijolo 


 Elês: 


Parede de uma vez e meia 


Parede de duas vezes 


O refechamento, tomado ou vincagem das juntas era indispensável nas 
paredes de tijolo porque evitava a infiltração da humidade, e proporcionava 
uma boa decoração. 


A solidez e a estabilidade estrutural de um edifício podia ser aumentada 
através da ligação das paredes extremas, sobre as quais assentava o viga- 
mento dos sobrados, por ferrolhos pregados numa das faces laterais ou na 
face superior das vigas. 


Referem-se ainda os pilares de tijolo que tinham habitualmente secção qua- 
drada, embora pudessem apresentar outras formas, obtidas com o auxílio 
de tijolos apropriados e cortados convenientemente. Mas, nestes casos, 
eram quase sempre rebocados e estucados para imitarem a cantaria, mais 
cara. [4] 


3.6 Paredes de adobe 


As paredes de adobe eram, em geral, construídas segundo as regras das 
de tijolo, com a particularidade de na sua grande maioria, apenas se empre- 
garem em construções pobres, fig. 2. 


De facto, estas paredes, correspondentes a uma solução ainda relativa- 
mente pouco conhecida, são constituídas por tijolos de barro cozido ao sol, 
ou em fornos, a temperaturas variáveis, baseados em matérias-primas muito 
diversas, onde a base de argila é composta com misturas arenosas ou cal- 
cáreas, que conferem ao barro características próprias [1]. 

Não podiam assentar directamento no solo devido ao risco de serem rapida- 
mente destruídas pela humidade ascendente, ao mesmo tempo que tornava 
as habitações insalubres. 

Habitualmente, os alicerces, que por esta razão tinham as suas faces supe- 
riores ligeiramente acima do terreno, eram de alvenaria corrente. 
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Fig. 2 Paredes de adobe — caso em Deixa-o-Resto (Santiago do Cacém) 


A partir deste nível, começava o assentamento dos adobes com argamassa; 
em geral, o tomado das juntas era feito de modo semelhante ao as paredes 
de tijolo. 


3.7 Paredes de taipa 


Como se refere na INTRODUÇÃO, a taipa corresponde a um sistema de 
construção milenar que perdurou entre nós até meados deste século na 
construção de edificações rurais. É formada por terra húmida comprimida 
entre taipais amovíveis de madeira, retirados depois de aquela estar com- 
pletamente seca, formando assim paredes ou muros homogéneos e monoli- 
ticos. 


A sua grande aplicação em construções rurais, principalmente no Algarve e 
Alentejo (mas também em algumas bolsas localizadas no Baixo Douro, na 
zona de Aveiro-Bairrada-Mira, na Estremadura (Vale do Tejo), e na raia 
junto a Castelo Branco [8]), deve-se ao facto de ser um material incombústi- 
vel, mau condutor de calor e, acima de tudo económico. 


No entanto, era facilmente atacada por roedores, e não era indicada para a 
construção de grandes edifícios, por razões de estabilidade face a acções 
sismicas. 


A taipa não podia assentar directamente sobre o solo porque as águas rapi- 
damente a desagregariam. Por isso, havia que fazer previamente a constru- 
ção de um pequeno alicerce ou sapata de alvenaria ordinária ou de tijolo 
com cerca de 0.25 m a 0.60 m de altura, acima do nível do solo, a partir da 
qual se dava início à construção da parede [4,6], fig. 3. 


Para se manterem os taipais na vertical, utilizavam-se caixilhos de 
desarmar, cada um deles constituído por um sarrafo de madeira horizontal e 
dois prumos, fixos em baixo por cunhas e em cima por uma corda mantida 
torcida com o auxilio de um pau, de forma idêntica ao cairo das serras de 
carpinteiro. Depois de colocados os taipais na sua posição, iniciava-se o 
fabrico da parede propriamente dito, colocando a terra ligeiramente húmida, 
em camadas de pequena espessura, de 10 cm em média, que se calcavam 
muito bem com maços apropriados, até que aquele valor fosse reduzido a 
metade. 

Após o enchimento da primeira série de taipais, estes eram desarmados, 
sendo colocados no seguimento da parede para que pudessem cobrir com- 
pletamente as faces laterais da parte inclinada do trabalho já concluído. 
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Fig. 3 Alicerces de uma parede de taipa — caso em Santiago do Cacém 


Assim que se completava a primeira fiada periférica da parede ou do muro, 
os taipais eram dispostos sobre a parte já construída de forma idêntica, 
tendo o cuidado de as juntas desta nova fiada não se sobreporem as da 


antecedente, apesar de serem inclinadas. 


Em geral, os vãos eram guarnecidos com tijolo, pedra ou madeira, sendo as 
vergas ou padieiras de madeira, na maioria das vezes, fig. 3. As paredes de 
taipa eram quase sempre rebocadas com argamassa de cal ordinária, tendo 
assim uma duração que podia atingir séculos: sendo para isso necessário 
que o reboco não fosse aplicado antes da taipa secar por completo (o que 
durava cerca de quatro meses) nem a parede devia ser rebocada com 
tempo húmido. A duração da taipa podia ainda aumentadar, se na confec- 
ção fosse humedecida com leite de cal. 

Após estar completamente seca, a taipa era finalmente revestida com arga- 
massa, conferindo-lhe uma grande duração. Mas, antes de se aplicar aquela 
camada de reboco, picava-se ligeiramente a sua superfície, molhando-a e 
chapando-a com cal e areia. 

Depois de rebocada, a parede podia finalmente ser calada. A caiação deve- 
ria ser repetida periodicamente, melhorando-se assim o estado de conserva- 
ção da parede. 


3.8 Paredes com armadura de palha 


Quando a argila era muito retractil, utilizava-se por vezes a palha como ele- 
mento aglutinador que, justamente, impedia aquele efeito, melhorando assim 
o comportamento higrotérmico da parede. 


A palha utilizava-se ainda noutras situações (ex. paredes de pedra seca) 
principalmente no norte do Pais, para vedar as juntas melhorando, assim, O 
seu comportamento térmico. 


3.9 Paredes de formigão 


Designa-se por formigão o aglomerado semelhante ao betão, em que a cal 
substitui o cimento. No QUADRO ll, apresenta-se a composição de alguns 
formigões, executados de forma semelhante à da taipa, principalmente no 
Alentejo. 


QUADRO Il - Composições de vários formigões, [4] 


SAIBRO | 
OU AREIA 
GROSSA 


ES 


(") Todos estes elementos eram amassados até se transformarem numa massa homogé- 
nia, que depois era colocada em moldes de madeira (taipais), semelhantes aos utilizados 
na confecção da taipa. 


3.10 Alvenaria de gesso argamassado 


Este tipo de alvenaria requeria cuidados especiais devido à rapidez com que 
se processava a presa, apenas podendo ser aplicada nos locais onde não 
fossem afectadas pela humidade do ar ou do solo. O gesso era amassado 
na altura da aplicação, nunca se interrompendo uma amassadura a meio, 
dado o perigo de endurecer rapidamente. 


As pedras eram primeiramente colocadas a seco no lugar que iam ocupar, 
sendo depois retiradas e colocadas junto à parede pela mesma ordem. 


Após a conclusão da amassadura do gesso, feita à mão e com colher de 
cobre, aplicava-se sobre a parede a argamassa necessária para se formar o 
leito e colocavam-se rapidamente as pedras nos lugares já “testados”. Com 
o gesso restante preenchiam-se os espaços vazios. 


3.11 Alvenaria de cimento romano 


Para se construir uma alvenaria de cimento romano havia que limpar muito 
bem os materiais e os locais de aplicação. No caso de alvenaria ser assente 
sobre muros, eles deveriam ser picados e muito bem varridos e lavados, 
para se garantir a necessária aderência. 


Exceptuando estes cuidados iniciais, o seu fabrico era relativamente sim- 
ples, sendo o cimento aplicado logo depois de amassado, dispondo-se 
de seguida as pedras e os tijolos sobre ele. Esta mistura era depois 
batida mas apenas durante os primeiros instantes, em que estava ainda 
fresca. 


A alvenaria assim construída, era conservada humida por meio de regas fre- 
quentes. 


3.12 Alvenaria de beton 


Definia-se como beton o aglomerado artificial que se obtinha pela mistura de 
saibro ou cascalho com argamassa hidráulica, geralmente de cimento Por- 
tland. 


Os principais impulsionadores das construções em betão foram os Roma- 
nos, existindo ainda hoje, em bom estado, algumas das suas obras tanto 
em Itália como em todas as regiões por eles dominadas. 

Mas, com a queda do Império Romano, desapareceu também a aplicação 
deste material, para voltar a aparecer em meados do século passado, princi- 
palmente na Alemanha e na Inglaterra. [6] 
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SUMÁRIO 


A construção do Aeroporto 
Internacional de Macau custou 
um pouco mais de um bilião de 
US& (cerca de 160 milhões de 
contos). O custo programado 
foi respeitado e a construção 
terminou dentro do prazo pre- 
visto. A presente nota descreve 
a maneira como foi preparado, 
organizado e conduzido todo o 
processo que permitiu a con- 
clusão do empreendimento 
com êxito. Termina com a ela- 
boração de considerações 
sobre alguns dos problemas 
mais dificeis que podem apare- 
cer em projectos desta natu- 
reza e apresenta as sugestões 
que com base na experiência 
colhida, se consideram uteis 
para O sucesso deste tipo de 
empreendimentos. 


SUMMARY 


Macau International Airport is a 
brand new airport built ofishore 
Macau at final cost of USS 
1 Billion. The whole project 
finished at cost and on time. 
This paper presents the work 
that has been done and how it 
has been prepared, program- 
med and conducted. The most 
difficult problems encountered 
during construction are high- 
lighted as well as the recom- 
mended short cuts to succeed 
in such a complex work. 


PREPARAÇÃO, CONDUÇÃO E 
ORGANIZAÇÃO DE UM 
GRANDE EMPREENDIMENTO 
DE ENGENHARIA 
(O CASO DO AEROPORTO 
INTERNACIONAL DE MACAU) 


António Diogo Pinto, Prof. IST 


1. INTRODUÇÃO 


Em 1988 estava concluído um Plano Director elaborado por um consórcio 
internacional liderado pelo Aeroporto de Frankfurt, o qual incluia uma 
empresa portuguesa, a Consulplano. 


A 18 de Janeiro de 1989 foi formada a CAM-Sociedade do Aeroporto Interna- 
cional de Macau, SARL, com o objectivo de construir e explorar o aeroporto. 


Em Abril de 1989 foi iniciado o projecto preliminar com vista à execução da 
obra e em Dezembro deste mesmo ano foram iniciados os primeiros traba- 
lhos. Estes trabalhos teriam um desenvolvimento lento e seriam mesmo inter- 
rompidos até Março de 1992, data em que retomaram o andamento normal. 


O aeroporto entrou em operacões em 9 de Novembro de 1995 e foi oficial- 
mente inaugurado em 8 de Dezembro de 1995 ou seja, seis anos após o pri- 
meiro início dos trabalhos e 44 meses após a retomada definitiva dos mesmos. 


Dado o elevado volume de trabalhos, a sua grande complexidade e a adver- 
sidade do meio natural onde os mesmos se desenvolveram foi necessário 
preparar uma metodologia pouco usual em situações análogas, mas que se 
revelou eficiente, conduzindo globalmente a resultados muito aceitáveis. Os 
sucessivos passos e fases da metodologia adoptada são apresentados nos 
pontos seguintes. 


2. REPARTIÇÃO DAS RESPONSABILIDADES ENTRE A CAM E OS 
ÓRGÃOS DE TUTELA 


A repartição de responsabilidades entre a CAM, concessionária para a cons- 
trução e exploração, e os Órgãos de Tutela, representados pela Autoridade 
de Aviação Civil de Macau (AACM), entretanto criada, foi regulada por inter- 
médio dum Contrato de Concessão. 


Nos termos deste Contrato de Concessão, a aprovação de todos os projectos 
do aeroporto, era da responsabilidade do Órgão de Tutela ou seja da AACM. 


A AACM dotou-se dos meios adequados para o efeito tendo recorrido a 
um excelente naipe de consultores qualificados, quer na área das matérias 
aeroportuárias quer na área dos projectos de índole puramente relacionada 
com a engenharia civil. 
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Fig. 1 Localização do Aeroporto 


3. ORGANIZAÇÃO E GESTÃO DOS PROJECTOS BÁSICOS 
DO EMPREENDIMENTO 


A CAM-Sociedade do Aeroporto Internacional de Macau, SARL, contratou 
um consórcio, denominado MIAC, liderado pela British Airporis, para a coor- 
denação da execução de todos os projectos do Aeroporto Internacional de 
Macau. Os projectos de natureza aeroportuária e alguns de engenharia civil 
foram elaborados pelo próprio consórcio e os restantes foram entregues, 
mediante ajuste directo, a consultores idoneos. 


Destes projectos ressaltam, pela sua importância e consequências, o pro- 
jecto da ilha artificial, adjudicado à empresa chinesa, “SCCCC - Shangai 
Communications Construction and Geaneral Contracting Corporation" e o 
projecto das fundações dos pavimentos da pista e dos caminhos de circula- 
ção que foi adjudicado à empresa portuguesa GRID. Este projecto foi poste- 
rrormente abandonado e substituido por outro, conducente a um método 
construtivo, mais simples, mais eficiente e mais económico. O esquema 
geral do empreendimento é apresentado na figura seguinte. 


4. CONTRATOS DE EMPREITADA 


4.1 Divisão dos projectos do empreendimento em grandes grupos para 
efeitos das empreitadas 


Para efeitos das empreitadas de construção, os projectos básicos do aero- 
porto foram divididos em três grandes grupos: 


a) a ilha artificial, pista e caminhos de circulação; 

b) os caminhos de ligação entre a ilha artificial e a área do terminal: 

c) o edifício do terminal, placa de estacionamento, torre de controlo, edifício de 
operações, equipamentos aeroportuários e equipamentos de aeronavegação. 


4.2 Escolha e contratação de empreiteiros 
Contrato n.º 1 - Ilha artificial, pista e caminhos de circulação 


A obra constituida pela ilha artificial, pista e caminhos de circulação na ilha 
artificial, foi objecto de contrato por ajuste directo, após uma pesquisa regio- 
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Fig. 2 Esquema geral do empreendimento 


nal sobre quais as empresas com capacidade, competência e meios para 
executar a obra e cumprir um programa muito exigente em termos de prazo, 
custo e qualidade. Foi seleccionado o Consórcio liderado pela “CHEC-China 
Harbour Engineering Corporation (RPC)” e constituído também pela “Zhen 
Wha Company (RPC)” e pela Tai Fok Wha, Ltd (Macau). 

Esta obra foi objecto de um contrato, por preço global fixo, na modalidade 
prevista na lei portuguesa, pelos valores seguintes: 

Ilha artificial 3038 milhões de patacas 
Fundações da pista e dos caminhos de circulação 1250 milhões de patacas 
Pavimento da pista e dos caminhos de circulação 250 milhões de patacas 


Nota: Uma Pataca = a vinte escudos 


Fig. 3 Vista geral da ilha artificial e dos caminhos de circulação 
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Contrato n.º 2 - Caminhos de ligação (pontes) entre a ilha artificial e a 
placa de estacionamento 


Esta obra constituída por duas pontes, uma de 1620 m de comprimento a 
Sul, e, outra de 700 m de comprimento a Norte, foi objecto de consulta a 
três consórcios tendo, após intensas negociações, o trabalho sido adjudi- 
cado a um consórcio liderado, pela “China Harbour Engineering Corpora- 
tion”, a que se juntaram a “Zhen Wha Company”, a Tai Fok Wha, Ltd e a 
empresa Construções Técnicas (Portugal). 

Esta obra foi também objecto de um contrato por preço global fixo no valor 
de 700 milhões de patacas, na modalidade dos contratos previstos na lei 
portuguesa. 


Fig. 4 Vista geral do caminho de circulação, em ponte, (lado Norte) e do topo Norte da ilha 
artificial 


Contrato n.º 3 — Edifício do terminal, placa, acessos, torre de controlo, 
equipamentos aeroportuários e equipamentos de aeronavegação 


Este contrato incluiu um conjunto de infra-estruturas de natureza essencial- 
mente aeroportuária que, na versão inicial, estava programado ser objecto 


Fig. 5 Vista geral do terminal (lado ar) e da placa de estacionamento 


Fig. 6 Vista geral dos Radares (principal e secundário) 


de oito contratos diferentes. A junção de todos estes contratos num só con- 
trato permitiu uma maior simplificação, racionalização e sobretudo uma clara 
alocação das responsabilidades durante a construção. 


A obra foi objecto de um concurso internacional e após prequalificação 
foram escolhidos três consórcios. Das negociações com os três consórcios 
saiu vencedor o Consórcio liderado pela empresa Soares da Costa, SA, O 
qual incluía também a Siemens “Plessey Systems, SA”. 


O valor do trabalho contratado foi de 1809 milhões de patacas e o contrato 
foi do tipo preço global fixo, na modalidade FIDIC (Federação das Associa- 
ções Internacionais de Empreiteiros e Construtores). 


A execução do contrato envolvia complexas operações de aprovisionamento, 
de construção e de ensaios pre-operação de índole essencialmente no domínio 
das disciplinas aeroportuárias. Por este facto decidiu-se afectar um “Programme 
and Construction Management” (PCM) exclusivamente a este contrato. 


5. COORDENAÇÃO E GESTÃO DO PROGRAMA, DO CUSTO E DA 
QUALIDADE 


5.1 Planeamento e gestão da construção 


Para o planeamento, fiscalização e controlo de custos e de prazos foram, 
após concurso público internacional, contratados dois consórcios, um para os 
trabalhos da área marítima, constituído pelas empresas “Scottwilson Kirkpa- 
trick”, Partex e “Pengest International” (PCM n.º 1), e outro, para os trabalhos 
da área terrestre, constituído pelas empresas americanas “Tams” e “Schall” e 
pela empresa portuguesa Profabril (PCM n.º 2). Os valores dos contratos 
foram, respectivamente, de 124 e 80 milhões de patacas, do tipo preço global 
fixo, e na modalidade prevista na lei portuguesa. Supletivamente, e em casos 
de dúvida, foi usada a filosofia FIDIC para a superação de potenciais conflitos. 


5.2 Programação da obra 


A programação de cada uma das três grandes parcelas da obra fazia parte 
integrante dos contratos estabelecidos com os empreiteiros (Contrato n.º 1, 
Contrato n.º 2 e Contrato n.º 3) e foi consolidada e harmonizada pelos gesto- 
res gerais da obra (PCMs). Para o efeito foi seleccionado e usado por ambos 
os PCMs, um programa de obra, actualizável todas as semanas (o programa 
“PRIMAVERA Project Planner" produzido por “PRIMAVERA SYSTEMS Ine.”). 


A obra concluiu-se no prazo programado de 44 meses. 
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5.3 Papel da CAM - Sociedade do Aeroporto Internacional de Macau, 
SARL 


A CAM, como dono da obra organizou-se no sentido de dar uma resposta 
eficiente às actividades da sua responsabilidade. Para o efeito criou uma 
estrutura aligeirada mas operacional. Dessa estrutura destacamos os três 
departamentos que tiveram um papel decisivo no êxito do empreendimento: 
a Comissão Executiva, responsável pela dinamização das acções e pelas 
decisões ao mais alto nível, o Departamento Administrativo e Financeiro, 
responsável pela mobilização, em tempo, dos recursos financeiros e pela 
elaboração, lançamento e acompanhamento dos concursos e subsequentes 
contratos, e, o Departamento de Engenharia, responsável pelo planeamento 
e controlo da execução da obra. 


- 
CAM's Engineering Dept | 


Project Engineer 


Progr. Const. Manager 
(Land Side) 


Progr. Const. Manager 
(Marine Side) 


Contract No. 1 | Contract No. 2 
| Artificial Island Conn. Tax. /Pav. 
Fig. 7 Organograma básico do Departamento de Engenharia da CAM 


6. FACTORES ADICIONAIS QUE CONTRIBUIRAM PARA O EXITO DO 
EMPREENDIMENTO 


A rapidez com que as obras foram executadas só foi possível graças a um 
conjunto favorável de factores, alguns dos quais já foram referidos nos pon- 
tos anteriores e outros, igualmente decisivos, que a seguir se enunciam. 
Destes destacamos, o processo de aprovação de projectos, o uso comum e 
quase generalizado da lingua inglesa e o uso dos mecanismos FIDIC, como 
instrumento contratual básico, nas empreitadas mais complexas e como 
meio supletivo, à lei portuguesa, para a superação de potenciais conflitos, 
nos contratos menos importantes. 


6.1 Simplificação do processo de aprovação dos projectos 


À Autoridade de Aviação Civil de Macau foi cometida a missão de aprovar 
todos os projectos, levados a cabo pela concessionária, na área da conces- 
são. Tal missão foi desempenhada com grande acerto e competência nos 
prazos requeridos pela programação para a maior parte dos projectos. Nos 
casos mais difíceis, em que tal não foi possível, a concessionária avançou 
sozinha para a execução dos mesmos, sob a sua exclusiva responsabili- 
dade. Foi o caso dos caminhos de circulação em ponte, das instalações 
para o combate a incêndios e do pavimento da pista e caminhos de circula- 


ção na ilha artificial, os quais só foram aprovados, muito tardiamente, em 
relação ao programa. Outros não foram mesmo aprovados, como foi o caso 
do pavimento. Nestes casos, o Governo remeteu para a concessionária a 
responsabilidade técnica e financeira correspondente aos riscos por esta 
assumidos, ao avançar com as respectivas obras sem autorização. 


6.2 Uso extensivo da língua inglesa 


Sendo um empreendimento conduzido por empreiteiros e consultores das 
mais variadas origens, culturas e línguas, necessário se tornou usar como 
língua veicular o inglês em todas as relações contratuais e de serviço. Nes- 
tas condições foi exigido o domínio da língua inglesa a todos os trabalhado- 
res com funções de chefia, ou que, desempenhavam funções de ligação 
entre as diferentes entidades empenhadas no empreendimento. 


6.3 Uso extensivo dos mecanismos FIDIC 


A variedade de empreiteiros e consultores de culturas e línguas diferentes 
aconselharam o uso de mecanismos contratuais e de relacionamento de 
natureza tão universal quanto possível. Os mecanismos previstos nos con- 
tratos FIDIC revelaram-se auxiliares eficientes e seguros, mesmo para, por 
similitude, arbitrar as lacunas dos contratos na modalidade prevista na lei 
portuguesa, como foi o caso dos contratos para as obras da ilha artificial e 
das pontes. 


7. CUSTO DO EMPREENDIMENTO 


O custo do empreendimento foi repartido nas seguintes rubricas: 


— Custos preparatórios; 

— Custos do Contrato n.º 1 ( Ilha artificial); 

— Custos do Contrato n.º 2 ( Caminhos de circulação em ponte); 

— Custos do Contrato n.º 3 ( Edifício do terminal, equipamentos aeroportuá- 
rios e equipamentos de aeronavegação); 

— Custos com a gestão da construção; 

— Custos com despesas gerais e outros. 


Os trabalhos a mais fixaram-se em 238 milhões de patacas o que repre- 
senta um pouco menos que 3% do custo total do empreendimento. 


No quadro seguinte são apresentados de modo resumido os valores corres- 
pondentes a esses custos: 


Custos dos trabalhos do Aeroporto e recursos disponibilizados 


Custos (milhões de patacas) Recursos (milhões de patacas) 
A — Custos preparatórios | - Capital social 
Obras 338 4000 
Projectos 306 
Subtotal B44 
B - Contrato n.º 1 — Ilha artificial J — Empréstimos bancários 
Valor base 4538 3200 
Trab. a mais 80 
Subtotal 4618 
C —- Contrato n.º 2 — Pontes K —- Ganhos de desenvolv. 
Valor base 700 2220 
Trab. a mais 20 


Subtotal 720 
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D - Contrato n.º 3 — Terminal 


Valor base 1809 
Trab. a mais 138 
Subtotal 1947 
E — Outros trabalhos 90 
F - Gestão da construção 204 
G - Despesas gerais 137 
H - TOTAL B360 L - TOTAL 9480 


Nota: uma Pataca = a vinte escudos (1 US$ = 8 Patacas). 


8. CONSIDERAÇÕES E RECOMENDAÇÕES FINAIS 


A análise dos pontos apresentados sugere que o êxito na execução do 
empreendimento se ficou a dever a um conjunto de opções tomadas desde 
o início pela Administração da CAM. 

Essas opções podem ser agrupadas, em função dos factores que as carac- 
terizam, do seguinte modo: 


* Em factores de simplificação, tais como: 
— simplificação e racionalização dos projectos; 
— simplificação do processo de aprovação dos projectos; 
* Em factores de selecção cuidada, tais como: 
— selecção dos empreiteiros com capacidade material e humana para 
executar os trabalhos; 
— selecção dos gestores do projecto e da construção (PCMs); 
* Em factores de normalização, tais como: 
— Uso extensivo da lingua inglesa em todas as comunicações de serviço; 
— uso da terminologia e do processamento FIDIC, mesmo nos casos em 
que os contratos base seguiam a lei portuguesa; 
* Em factores de organização, tais como: 
— a adopção de uma organização funcional muito simples; 
— a adopção de uma instrumentação simples, eficaz e com provas dadas 
para o controlo do progresso da obra (o “PRIMAVERA project planner”). 


A adopção das medidas acima referidas criou uma base de entendimento e 
confiança entre todos os intervenientes do processo. Por outro lado todos 
conheciam perfeitamente o seu papel específico e as interfaces com tercei- 
ros. Finalmente todos se sentiam continuamente controlados e o respectivo 
trabalho avaliado e recompensado ou penalizado. 

Estas foram as bases para o sucesso na execução da obra de engenharia 
acima apresentada. 

Espera-se que esta nota seja útil para os engenheiros que tiverem que se 
dedicar a trabalhos desta natureza. 


BIBLIOGRAFIA 


Aeroporto Internacional de Macau (1989-1995). Projectos e Relatórios de 
Obra, Documentos de trabalho não publicados. Macau. 


International Civil Aviation Organization (1990) “International Standards and 
recomended practices: Aerodromes — Annex 14 to the Convention on Inter- 
national Civil Aviation”, Vol. |, Aerodrome design and operations. 


Jan-96 


RESUMO 


Neste artigo  analisam-se 
aspectos relevantes de termo- 
dinâmica e cinética da corrosão 
pelos gases e sais fundidos. 


ABSTRACT 


Thermodynamic and Kinetic 
features of high temperature 
corrosion associated with gases 
and moltem salt contaminants 
are toriefly described. 


CORROSÃO A ALTA 
TEMPERATURA 


César A.C. Sequeira 


1. INTRODUÇÃO 


O estudo da corrosão a alta temperatura constitui um tema de grande 
importância científica e tecnológica, como o atestam os inúmeros artigos e 
livros publicados sobre este assunto, bem como as diversas reuniões inter- 
nacionais a ele consagradas [1-8]. Essa importância resulta, em grande 
parte, do emprego crescente de materiais de engenharia em equipamentos 
que operam a alta temperatura, aos quais se exige boa resistência mecã- 
nica e anti-corrosiva a temperaturas cada vez mais elevadas. De entre os 
equipamentos sujeitos a esta forma de corrosão podem citar-se reactores 
nucleares, geradores magnetohidrodinâmicos, caldeiras de vapor, permuta- 
dores, turbinas de gás e de vapor, motores diesel, pré-aquecedores de ar, 
economizadores, certas peças de fornos, etc. 


Para além do oxigênio, que é o principal agente corrosivo nos meios reacti- 
vos a alta temperatura, são também de destacar os seguintes agentes 
agressivos: 


— Enxofre e gases contendo enxofre (importante na indústria petrolífera e 
petroquímica, por exemplo); 
— Hidrogênio (indústria química, descarbonetação de aços, etc.); 


— Halogéneos e compostos halogenados (particularmente agressivos em 
atmosferas muito quentes); 


— Vapor de água (centrais eléctricas, sistemas aquosos sob pressão, etc.); 
— Azoto e amónia (origina endurecimento superficial dos aços); 


— Substâncias fundidas (processos electrometalúrgicos, energia nuclear, 
engenharia quimica, etc.); 


— Carbono e gases contendo carbono (reactores nucleares, etc.); 
— Cinzas (motores diesel, caldeiras, etc.); etc. 


À corrosão nestes sistemas reveste-se de aspectos muito diversos sendo 
extremamente difícil, senão mesmo impossível, classificar e definir em ter- 
mos teóricos as várias observações e fenómenos. Há, no entanto, dois tipos 
de sistemas que se encontram com mais frequência e que têm sido estuda- 
dos com mais profundidade: — trata-se dos sistemas do tipo metal/gás e do 
tipo metal/sal fundido. Há também, como sucede com qualquer outro pro- 
blema químico, dois aspectos básicos a considerar nestes tipos de siste- 
mas: a termodinâmica dos processos e a cinética dos processos. Neste 
artigo, serão discutidos alguns aspectos termodinâmicos e cinéticos da cor- 
rosão pelos gases e sais fundidos. 


Os sistemas do tipo material/metal fundido têm vindo a ganhar importância 
em virtude do uso crescente de metais fundidos como refrigerantes líquidos 
em reactores nucleares [9, 10] e como componentes de baterias [11] e dou- 
tros sistemas de alta temperatura utilizados em estudos de termodinâmica e 
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cinética de ligas metálicas fundidas e outros líquidos iónicos [12], com vista 
à sua adequada aplicação em processos de refinação industrial, entre 
outros [13, 14]. Apesar da presente relevância destes sistemas, falta de 
espaço não nos permite fazer um estudo, mesmo que resumido, dos aspec- 
tos mais relevantes que lhes dizem respeito. 


2. OXIDAÇÃO A ALTA TEMPERATURA 


2.1 Reacções metal-gás 


A maior parte dos metais encontram-se no meio ambiente associados a 
não-metais; isto significa que o seu estado mais estável é o que resulta da 
transferência de electrões do metal para o não-metal de forma a constituir- 
se o composto. O termo “oxidação metálica” aplica-se a qualquer processo 
que prive o metal de electrões quando ele é exposto a um meio ambiente 
“seco”. O crescimento de filmes na superfície metálica, é um aspecto essen- 
cial da oxidação, o qual é influenciado por características termodinâmicas, 
químicas, físicas e mecânicas do filme. Esse filme actua como uma camada 
barreira entre o metal e o oxidante, e o prosseguimento do seu crescimento 
assim como a reacção metal-oxidante é função da forma como os catiões 
metálicos, electrões e aniões do oxidante são transportados através do 
filme. Isto ilustra-se esquematicamente na Fig. 1. Nela se pode observar 
que o filme de óxido serve simultaneamente como um condutor iônico (elec- 
trólito), um condutor electrónico, um eléctrodo no qual se reduz o oxigénio, e 
uma barreira de difusão através da qual iões e electrões devem migrar. 
Efectivamente, a oxidação metálica M+ 140,» MO consiste de dois pro- 
cessos parciais complementares, a oxidação de M a M* na interface metal- 
óxido e a redução simultânea de O, a O” na interface óxido-gás, e claro que 
só se processará se electrões e iões forem conduzidos através do filme. 


M— A?" + 2e- 
ce-lis 1205 + 2e- — 02 


M + 20» —— MO 


Oxigênio 


Fig. 1 Processos electroquímicos que ocorrem na oxidação gasosa. 


Na prática corrente, os metais são geralmente impuros, contendo impurezas 
metálicas e O,H,N,C,S, etc. As atmosferas gasosas que lhes estão associa- 
das são também, em geral, muito complexas, contendo O, N.. H,, CO,, H,O, 
etc. Por conseguinte, os processos de oxidação que ocorrem nos sistemas 
do tipo metal-meio gasoso a alta temperatura são processos químicos extre- 
mamente complicados. 

Na Fig. 2 ilustra-se a grande variedade de fenómenos químicos e físicos 
que se processam na rede metálica e à superfície dum metal num meio 
gasoso típico a alta temperatura. Os vários fenómenos sistematizam-se em 
cinco zonas de actividade: metal, interface metal-óxido, óxido, interface 
óxido-gás, e gás. 


As reacções químicas que se podem verificar num sistema do tipo metal- 
gás a alta temperatura podem classificar-se em nove tipos como se indica e 


exemplifica a seguir: 


1. Oxidação directa, nitretação ou sulfuretação; no caso da oxidação directa 


ter-se-á por exemplo: 


Zr+ 0, = Zro0, 


Metal, M Óxido, MO; | | Das 


(1) 


Q 1 

E 

ais 

ada 

PAES 

lona À Metal lona É Óxido | lona E Cas 

Solução ou evolução de Difusão de ONHMLC eS| Difusdo qososo de O, Na, 
MH SE. Nuciecção e crescimento Ha, 50, etc. para o su- 
Formação de CH e | do gxido. | perficie. 
Precipifoção de oxidos, Reocçoo de M,5S eC com 
mitretos e hidretos O. 


Migração electrónico. 
Redução de óxidos delta 
valêncio pelo M, 
Reacções interóxidos 


Zona B Interface Metal -Óxido lona D Interface Óxido -Gás 


Tronsferêncio electrónico Adscrçdo de O. NH H50. 
Solução de MO NHCEC,S Solução de ON HS ec. 

nos óxidos ou em M Desorção de CO, SOs, ete. 
Reacção de M com O Formeção de âxdo 

Formoçoo de oxianives Nucieoção e crescimento de dxrdo 
Reocçoo ce C e S comô, Dissociação do Gxido. 


Volatilizoção de Me MO, 
Formoção de locunas 
Tronsterêncio electrônica 


Fig. 2 Actividade quimica de um metal numa atmosfera activa a alta temperatura, (reac- 


ções tipo óxido). 
2. Formação directa de óxidos voláteis: 


E ? 
Mo+-—0, = “(MoO,) 
) a 3 1 


3. Dissolução directa de O.N.H ou S no metal: 


Zr+0, =(Zr)j0, (solução sólida de O.) 


4. Reacção com H,O: 


Zr+2H,0=Zr0O,+2H, (solução sólida e/ou gás) 


5. Oxidação directa de óxidos de valência baixa: 


Ea mê 
2Fe,0,+-0,=3Fe,0, 


(2) 


(3) 


(4) 
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6. Formação de óxidos complexos: 


Fe0+ Fe,0, = Fe.O, (6) 


7. Redução de óxidos de valência elevada pelo metal: 


Fe,0,+ Fe=4Feu (7) 


8. Redução de óxidos pelo C ou 5: 
Fe0O+C = Fe+CO (8) 


9. Volatilidade dos óxidos: 


O) =| HF 

3Wo, =(WO,), (9) 
A ocorrência de uma determinada reacção depende da temperatura, pres- 
são gasosa, caudal gasoso e propriedades físicas e químicas do sistema 
metal-óxido. Observe-se que nos exemplos dados acima se evitaram deta- 
lhes concernentes às diversas reacções, é por isso que não se incluiram as 
reacções que são acompanhadas pela formação de lacunas. 
Num sistema do tipo liga metálica binária-gás as reacções que podem ocorrer 
são ainda mais numerosas e complexas, como é óbvio; resultam então com- 
postos binários intermédios (12CaO-7ALO,), soluções sólidas com miscibili- 
dade total (FeO-NiO), soluções sólidas com miscibilidade parcial (ALO,-Fe,O,), 
óxidos duplos do tipo espinela (ZnALO,), etc. A título de exemplo, vejamos as 
reacções que é possível observar aquando da oxidação ao ar a alta tempera- 
tura, das ligas de ferro-crômio: 


1. Oxidação dos elementos da liga 
Fe+0, — FeO — Fe 0, + Fe,0, (1 0) 


2. Reacções no estado sólido 
(a) Combinação entre óxidos 


Fe0O+CrO, > FeCro, (12) 
FeCrO, + Fe,0, > Fe(Fe,Cr),0, (13) 
Fe,0, + CrO, > (Fe.Cr),0, (14) 
(b) Difusão do oxigénio 
FeO+ 0 — FeO, > Fe,0, (15) 


2Fe(Fe.Cr),0,+ a O, > K(Fe.Cr),0, 


(16) 
(c) Difusão do metal 
4(Fe.Cr)O,+(Fe.Cr) > 3Fe(Fe,Cr),0, 
Fe,0, + Dear Preo (17) 


FeoO, 4+Cr— amam 2Fe +CRU, (18) 


Fell, = Foi) 


Está claro que no caso da liga ser constituída por múltiplos elementos e, 
além disso, o meio gasoso constar de uma mistura de gases, é possível 
arranjar diversas classificações, as quais englobarão inúmeros tipos de 
reacções. 


Sob o ponto de vista do corrosionista, são os sistemas em que ocorrem 
reacções extensas entre metais e gases activos a temperaturas elevadas, 
que tem mais interesse considerar, pelo que esta exposição se confinará a 
eles. 


Este campo da metalurgia química é particularmente interessante pelo facto 
da actividade química de um metal nas várias atmosferas agressivas ser um 
factor importante na extracção, processamento e utilização do metal. 
A obtenção de melhores metais e ligas metálicas para usos a alta tempera- 
tura em atmosferas reactivas especiais, e de metais com propriedades fisi- 
cas determinadas, tém sido um dos principais objectivos dos metalurgistas 
nas ultimas décadas. As recentes indústrias baseadas na energia nuclear, 
exploração espacial e electrónica molecular apresentam muitos problemas 
de materiais novos. 


Os cientistas espaciais pretendem materiais com elevada resistência mecá- 
nica que aguentem temperaturas acima de 2600ºC em meios gasosos oxi- 
dantes em regime turbulento. Os engenheiros nucleares desejam materiais 
com propriedades nucleares adequadas que resistam a atmosferas aquosas 
e de vapor a temperaturas da ordem dos 800ºC e a altas pressões durante 
longos anos. Outros procuram materiais que resistam a atmosferas de dió- 
xido de carbono a 800ºC, também a altas pressões durante longos anos. Há 
também os técnicos de corrosão cujo objectivo é proteger os automóveis, 
edificios e outras estruturas metálicas contra a corrosão em atmosferas 
industriais e marinhas. 


Os pioneiros da oxidação pensavam que o mecanismo básico das reacções 
metal-gás consistia numa difusão simples de iões ou átomos do metal e do 
oxigénio através de um filme ou película uniforme de óxido recobrindo a 
superficie metálica. Hoje sabe-se que o mecanismo é muitíssimo complexo. 
Para além da difusão, muitos outros factores determinam a natureza e 
extensão das reacções metal-gás. Alguns desses factores são: — a estrutura 
metálica, estrutura do óxido, difusão metálica, difusão do oxigénio ou do 
metal nos limites de grão ou à superfície, volatilidade do metal e do óxido, 
solução de oxigénio no metal, empolamento, fissuração, quebra em lascas, 
sinterização, nucleação do óxido, condutância electrónica, adesão do óxido 
e plasticidade do óxido. 

As reacções metal-gás podem auxiliar-nos na interpretação da corrosão em 
fase líquida, domínio que é ainda mais complexo. Com essa ajuda, é possi- 
vel identificar prontamente as zonas metálicas onde a corrosão localizada se 
instaura usando-se métodos electrónicos ópticos no estudo topológico da 
superfície do metal. 


De entre estes métodos destaca-se a microscopia electrónica de varrimento 
(SEM), a microscopia electrónica de transmissão (TEM), a sonda electrónica 
de microanálise (EPMA), a espectroscopia de fotoelectrões X (XPS), a 
espectroscopia de absorção de raios X de estrutura fina (EXAFS), a espec- 
troscopia de fotoelectrões de ultravioleta (UPS), a espectroscopia de elec- 
trões Auger (AES), a espectroscopia transformada de Fourier de |.V. (FTIR), 
a difracção de electrões de baixa energia (LEED), a espectroscopia de 
massa de iões secundários (Static SIMS), a espectroscopia de retrodifusão 
de Rutherford (RBS), a espectroscopia de Raman (RAMAN), a espectrosco- 
pia de raios X induzida por protões (PIXE) e a espectroscopia de perda de 
energia de electrões lentos (HREELS) [15-18]. 
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2.2 Termodinâmica da oxidação 


Uma das primeiras questões com que o corrosionista depara consiste no 
estabelecimento das condições termodinâmicas das reacções. A energia 
livre de formação dos óxidos, nitretos, carbonetos, sulfuretos e hidretos usa- 
se no solucionamento desta questão. Como é de prever, a energia livre de 
formação destes compostos varia muito. Consequentemente podemos recor- 
rer às energias livres de formação para classificar as condições de formação 
de óxidos, sulfuretos, hidretos, etc. 

Como se sabe, a energia livre de formação de um composto AGº, fornece 
uma medida da afinidade química a temperatura e pressão constantes. Para 
precisar isto, considere-se a formação do óxido de niquel a partir do níquel 
e duma mole de oxigênio: 


NH O, = 2NºO (19) 


Designemos por P(O,) a pressão do oxigénio no equilíbrio, ou tensão de dis- 
sociação do óxido a dada temperatura. Neste caso, a constante de equilibrio 
é dada por K=I/P(O.); uma vez que AG'º e K estão ligados pela relação 
AG" =-RT InK = «4,575 logk, pode, pois escrever-se: 
AGº=4,575 T log P(O,) (20) 
Nestas equações as pressões exprimem-se em atmosferas e a variação de 
energia livre em calorias. Ora, a quantidade 4,575T log P(O,) é igual à dife- 
rença entre o potencial químico de uma mole de oxigénio combinado no 
estado de óxido na presença do metal que o originou, à temperatura T, e o 
potencial químico de uma mole de oxigénio gasoso, a 1 atmosfera e à 
mesma temperatura. Por conseguinte, AGº caracteriza o potencial quimico 
do oxigénio no óxido e -AGº a afinidade padrão do metal para o oxigénio. 
Quanto mais negativo é o valor de AG”, mais favorável é a formação do 
Óxido, e mais básico é o metal. 
Encontram-se na literatura numerosos dados sobre a energia livre de forma- 
ção de diferentes óxidos, sulfuretos, hidretos, etc, [19-24]. Mas é mais fácil 
usar diagramas como o da Fig. 3 relativo aos sulfuretos [19], no qual a tem- 
peratura é colocada em abcissas e a energia livre de formação padrão por 
mole de enxofre é colocada em ordenadas. Cada linha do diagrama é a 
expressão gráfica da equação AG' = AHº-TAS”, considerando AH" e AS 
independentes da temperatura. A equação de AG” em função de T é por- 
tanto uma recta cuja ordenada na origem é o calor da reacção e que tem 
por declive a variação de entropia. Quando AH" e ASº variam com a tempe- 
ratura, os termos introduzidos pelo calor específico compensam-se e a 
equação AG' = f(T) permanece praticamente uma recta. No entanto, quando 
ha uma mudança de estado, esta é acompanhada por uma alteração no 
valor AS”. Assim uma fusão ou volatilização do sulfureto diminui o declive da 
recta e uma fusão ou volatilização do metal aumenta-o. 


Estes diagramas, conhecidos por diagramas de Ellingham ou de Ellingham- 
Richardson, apresentam um grande interesse. Por um lado, a temperatura 
elevada T A Sº cresce e por conseguinte a afinidade dos metais para o oxi- 
gênio diminui com o aumento da temperatura; aumenta a tendência dos sul- 
furetos, óxidos, etc. a decomporem-se. A intersecção das rectas com o eixo 
das abcissas dá a temperatura para a qual a tensão de dissociação dos óxi- 
dos, sulfuretos, etc., é igual a uma atmosfera. Por outro lado, estes diagra- 
mas são equivalentes a conjuntos de equações com dois parâmetros. Tem- 
se, com efeito, AG" = A-BT e log K = C/T+D. Estas equações são preferíveis 
às que, incorporando os coeficientes de todos os calores específicos, dão 
lugar a cálculos complicados e a uma precisão ilusória. 


Quando se pretende comparar a aptidão para a oxidação de diferentes 
metais é possível, em primeira aproximação, recorrer aos valores AHº em 
vez de AGº. É que as variações de AHº pouco diferem das de AGº quando 
se passa dum metal para o outro. Uma escala aproximada das afinidades 
dos metais para o oxigénio pode pois obter-se classificando os metais de 
acordo com os valores crescentes dos calores de formação dos seus óxi- 
dos. A facilidade de formação do óxido, por conseguinte, depende da afini- 
dade do metal para o oxidante, da temperatura da reacção e da pressão (e 
actividades) do sistema metal-oxidante. 


A aplicação de valores termodinâmicos no cálculo das condições de viabili- 
dade das reacções de oxidação deve fazer-se com certas reservas. Muitas 
reacções envolvem (a) a formação de gases que são romovidos da zona de 
reacção ou (b) a formação de gases que condensam aquando do arrefeci- 
mento e são removidos da zona de reacção. Nestes casos as reacções 
podem verificar-se mesmo em condições termodinâmicas desfavoráveis por- 
que o equilíbrio nunca chega a ser atingido. Um exemplo interessante do 
que se acaba de afirmar é a reacção que se desenrola entre o carvão e o 
hidrogénio nos reactores de fluxo a altas temperaturas, caso em que todos 
os hidrocarbonetos são instáveis. 
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Fig. 3 Diagrama de energias livres de sulfuretos metálicos. 
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Apesar das limitações da termodinâmica, é um facto generalizado que, em 
todas as questões referentes ao ataque de metais ou ligas metálicas por um 
ou mais oxidantes, é costume abordar essas questões com uma análise ter- 
modinâmica tão completa quanto possível das condições necessárias para 
que o ataque seja viável. Evitam-se assim muitos estudos cinéticos que, 
além de serem complexos, se tornariam muito morosos [25]. 


2.3 Cinética da oxidação 


Quando um material metálico é sujeito a uma atmosfera oxidante, há forma- 
ção de um filme que vai ditar, de acordo com as suas características, a pos- 
sibilidade do processo de oxidação prosseguir com uma velocidade que 
será fortemente dependente dos fenómenos de transporte através do refe- 
rido filme. 

Um dos principais objectivos dos estudos das reacções de oxidação é o 
estabelecimento de leis de velocidade empíricas e, sempre que possivel, 
leis de velocidade teóricas baseadas nos modelos físicos explicativos da 
cinética das reacções. As leis de velocidade empíricas permitem extrapolar 
velocidades de reacções em novas condições experimentais enquanto que 
as leis de velocidade teóricas conduzem, em adição, ao mecanismo da 
reacção. O conhecimento completo do mecanismo da reacção permite o 
estabelecimento de predições acerca doutras reacções. 


A teoria com mais aplicação às reacções sólido-gas a alta temperatura, 
especialmente nos casos em que as camadas de oxidação resultantes são 
muito espessas, é a de Wagner [26]. Wagner admite que o transporte de 
matéria através da camada de óxido se efectua sob a forma de iões e de 
electrões, pela acção dum gradiente de concentração devido a uma varla- 
ção de composição do óxido entre as interfaces metal/óxido e óxido/gás, 
e/ou por influência dum gradiente de potencial eléctrico devido a diferenças 
de concentração de carga, em particular nas interfaces. Por conseguinte, 
segundo Wagner é a difusão da rede propriamente dita que é a etapa deter- 
minante do processo de oxidação. 


Tendo em atenção os estudos acabados de citar e outros discutidos por 
Kofstad [27], Bénard [28], Hauffe [29], etc., é possível classificar os vários 
processos controlantes da velocidade das reacções de oxidação em quatro 
grupos principais: (i) difusão do gás da reacção ou dos produtos da reacção 
nos sentidos gás/superfície metálica e superfície metálica/gás, (ii) reacções 
nas interfaces metal/gás ou óxido/gás, (ili) processos de difusão do metal, 
gás ou electrões através do óxido ou na fase metálica e (iv) processos fisi- 
cos. Para estes grupos principais de processos foi possível estabelecer leis 
de velocidade teóricas. Destas destacam-se os cinco tipos principais seguin- 
tes: — (Ii) linear, (il) parabólica, (Iii) cubica, (iv) logaritmica directa e (v) loga- 
rítmica inversa (Fig. 4). Em casos extremos a cinética da reacção de oxida- 
ção pode alterar-se e passar a observar-se uma lei rectilinea. A lei 
logaritmica inversa, que constitui o outro caso limite, verifica-se a temperatu- 
ras baixas às quais o movimento iónico só é possível na presença de um 
gradiente de potencial eléctrico. Outras leis logaritmicas menos gerais inter- 
pretam-se como devidas a formação de imperfeições mecânicas nos filmes 
de oxidação tais como bolhas, poros e/ou fendas. Uma lei cúbica indica uma 
camada protectora. 


As leis de formação de filmes de oxidação observadas a altas temperaturas 
são a linear, a parabólica e a paralinear ou parabólica mista. Contudo, diver- 
sas complicações inerentes à oxidação dos metais e ligas metálicas introdu- 
zem factores desconhecidos que resultam em desvios consideráveis das leis 
simples de crescimento de óxidos aplicáveis à formação de filmes coeren- 
tes. 


2.3.1 Lei de crescimento linear 


Se o peso W de produto de oxidação formado por unidade de superficie do 
metal base é directamente proporcional ao tempo de oxidação t, isto é se 
W=K,, então a representação gráfica de W em função de t traduz-se por 
uma linha recta. O crescimento linear ocorre desde que a formação dos pro- 
dutos de corrosão seja determinada pela velocidade de provimento de oxi- 
dante à superfície da amostra metálica, por exemplo a baixas pressões de 
oxigénio. Isto significa que a disponibilidade de iões metálicos não é o factor 
determinante da velocidade, sendo esta efectivamente controlada pela exis- 
tência de oxidante necessário à formação do filme de óxido. 
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Fig. 4 Leis da velocidade de oxidação geralmente observadas. 


A ionização do oxidante, a quimisorção e/ou a nucleação superficial são 
outros factores controlantes do crescimento linear. Hauffe e Rahmel [30] 
observaram que a oxidação do Ni pelo vapor de enxofre na zona de pres- 
sões 0,01-0,5 mmHg segue uma lei linear, no entanto, à pressão de 
1 atmosfera a lei de crescimento é parabólica [31]. 


A oxidação linear pode ocorrer por outras razões. Um óxido descontínuo, 
com fissuras, poros, ou um cuja formação é acompanhada pela introdução 
de grandes tensões compressivas no óxido (originando fendas e subcama- 
das não aderentes e pulverulentas) pode expôr o metal continuamente ao 
gás atacante a altas temperaturas ou então o óxido pode romper-se ou des- 
prender-se aos bocados (em pó) à medida que se forma. Esta rotura espon- 
tânea ocorre sempre que a espessura dum óxido sob tensão mecânica 
ultrapassa um certo valor critico. A velocidade de oxidação passa então a 
ser controlada pela camada inicial sob tensão de espessura constante, 
resultando duma sucessão de velocidades parabólicas cujo valor médio jus- 
tifica o comportamento linear do processo global. 


2.3.2 Lei de crescimento parabólico 


O crecimento de filmes de acordo com a relação W* = K, t + const. aplica-se 
a muitas reacções de formação e crescimento de camadas oxidantes a alta 
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temperatura. A difusão através da rede por via dos seus iões e electrões é o 
factor controlante da velocidade. Claro que no caso de camadas muito 
espessas o gradiente de potencial eléctrico é desprezável, apenas intervindo 
o gradiente de concentração, isto é, é a difusão iónica que controla a veloci- 
dade. 

Uma interpretação simples da lei parabólica é a seguinte: — Comparem-se 
as duas situações representadas na Fig. 5. Os limites do gradiente de con- 
centração através do óxido têm os mesmos valores em ambos os casos, de 
acordo com a regra das fases (que determina que a composição de um 
composto seja fixa desde que a temperatura, e a composição e pressão do 
gás circundante o sejam). Por conseguinte, o tempo necessário para um ião 
metálico difundir do metal para o gás no caso Il será quatro vezes maior 
que no caso |; não só o gradiente de potencial químico se reduziu a metade 
do valor que tinha em I[Au(II/AX=Ay(I/2AX| como o ião tem também um 
trajecto duplo a percorrer [espessura do filme em II (2X) é dupla da espes- 
sura do filme em I (X)] . O grande mérito desta “derivação” é que nos da 
uma imagem física da lei de crescimento parabólico da difusão. Wagner [32] 
deduziu uma equação quantitativa da oxidação parabólica, conhecida por 
Equação de Wagner: 


Ii 
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Fig. 5 Crescimento de filmes segundo a lei parabolica. 
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O termo da direita chama-se constante racional de velocidade, Kr (em equiv. 
gr/'em seg.), de modo que dM/Adt = K/X, onde X é a espessura de óxido 
instantânea (cm) no tempo t(seg) e M é o número de equivalentes de óxido 
formados no tempo t; (M**) é o número de moles de catiões por cm”, z+ e z- 
são as cargas dos catiões e aniões do filme, D e D são as autodifusivi- 
dades dos catiões e aniões, (M), (MM) e (MG) são as actividades do metal 
no óxido, na interface metal/óxido e na interface óxido/gás, respectivamente. 
Claro que X=const M, visto que M é o número de equivalentes de óxido ins- 
tantâneos no tempo t. Portanto, dM/dt = K, A/X = consUM, de modo que inte- 
grando vem Mº = const. tou Wº = K1 (lei de crescimento parabólico, com W 
em mg de óxido formado por cmí da área A do meial). K, designa-se por 
constante de velocidade parabólica e claro que K,=const. K,. Pode escrever- 
se K,=2xI0º dE K, onde d é a densidade do óxido (g/em'*), E é o peso do 
seu equivalente grama (9) e K,, se exprime em mg'/cm' seg. K, pode também 
derivar-se em função dos números de transporte e das mobilidades dos 
catiões e aniões da camada de óxido, e foi calculado para diversos óxidos, 
sulfuretos, etc., para os quais K, se obteve experimentalmente, verificando- 
se plena concordância entre os valores calculados e os teóricos. 


À variação da difusividade com a concentração ou da velocidade de cresci- 
mento de camadas com a mesma espessura a temperaturas diferentes 
devida a modificações nos seus limites de composição não invalida a lei 
parabólica observada em condições particulares. Contudo, a variação da 
composição da fase metálica durante a oxidação (por exemplo por redução 
ou eliminação de um componente da liga) invalida as condições limites de 
validade da equação de crescimento parabólico. 


2.4 Mecanismos da oxidação 


O estudo cinético das reacções metal-gás é uma forma prática de avaliar a 
reactividade química dos metais. Por seu lado, o estabelecimento das leis 
de velocidade teóricas conduzem, em princípio, ao mecanismo das reac- 
ções. 

As três teorias principais interpretativas da oxidação dos metais são as de 
Cabrera e Mott [33], Wagner [26] e Gulbransen [34]. Cabrera e Mott suge- 
rem que o factor controlante do crescimento de filmes delgados é a existên- 
cia de um campo eléctrico poderoso (cerca de 10” volt/cm) através do filme. 
Outros factores tais como a quimisorção [35], a transferência de electrões 
na interface [36] e a difusão associada à formação de cavidades [37] foram 
aventados. Kubaschewski [38] sugere que o modelo de Evans [37] é o que 
mais provavelmente se aplica a filmes delgados. 


Os estudos da estrutura cristalina dos metais feitos antes, durante e depois 
das reacções são um auxiliar poderoso na descrição e classificação das 
reacções de oxidação. Os pioneiros destes estudos bem como dos de 
difracção electrónica pensavam que essa ferramenta por si só era suficiente 
para a compreensão dos processos corrosivos. Hoje sabe-se que a análise 
completa da estrutura cristalina de um metal e dos seus filmes superficiais é 
apenas uma parte, embora importante, dum problema extremamente com- 
plexo. 


Para a elucidação dos mecanismos da oxidação, muitos mais aspectos con- 
vém considerar, para além dos fenómenos de estrutura cristalina associados 
aos processos de oxidação. Em particular, os seguintes parâmetros devem 
ser tidos em consideração [25]: 

. Função de trabalho electrónico no metal 

. Energia de remoção dos átomos metálicos 

. Solubilidade do oxigénio, enxofre ou hidrogénio no metal 

. Volatilidade do óxido e metal 

- Temperaturas de fusão do óxido e metal 

. Formação de óxidos multiplos 

. Transformação da estrutura cristalina do óxido e metal 

- Coeficiente de expansão do óxido relativamente ao metal 

. Condução iónica e electrónica dos óxidos. 


A condução iónica e electrónica dos óxidos é uma propriedade que desem- 
penha um papel crucial nos processos de oxidação. Os óxidos metálicos 
(sulfuretos etc.) são usualmente cristais iónicos imperfeitos, pois caso fos- 
sem perfeitos funcionariam como isolantes eléctricos impedindo o cresci- 
mento dos filmes de oxidação, logo após o seu aparecimento na superfície 
dos metais. Esse crescimento, característico das reacções sólido-gás a alta 
temperatura, ocorre, por conseguinte devido à natureza defeituosa desses 
filmes e à consequente difusão que se verifica através deles. O estudo dos 
defeitos ou imperfeições pontuais das redes cristalinas foi iniciado cerca de 
1930 por Schóttky, Wagner, Jost e Frenkel. Estes cientistas fizeram um 
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estudo geral do estado cristalino. Segundo eles, existe num cristal qualquer, 
em equilíbrio termodinâmico, um certo número de lugares da rede vazios a 
que chamaremos lacunas ou, por outras palavras, um certo número de pon- 
tos da rede em que o átomo ou ião correspondente está ausente. Schóôttky 
e Frenkel interpretaram sobretudo a estrutura defeituosa dos compostos 
iónicos estequiométricos (ex. NaCl, CaF,. AgBr e PbCI,) [39], enquanto que 
Wagner estudou a não-estequiometria dos semicondutores. Dentro do tema 
que estamos a tratar interessa especialmente abordar o caso dos óxidos 
não estequiométricos e dos seus defeitos, assim como as consequências 
destes no processo de difusão, pelo que lhe dedicaremos a subsecção 
seguinte. 


2.4.1. Condução nos compostos não estequiométricos 


O mecanismo da condução nos compostos não estequiométricos é um caso 
de estudo particularmente interessante para compreender a condução nos 
óxidos a altas temperaturas, porque a alta temperatura, poucos óxidos são 
estequiométricos. A não estequiometria é de facto um dos mecanismos pelo 
qual aparece a semiconductividade. À não estequiometria aparecem asso- 
ciados defeitos pontuais, logo este fenómeno aumenta simultaneamente as 
conductividades electrónica e iônica [40, 41]. 


Os desvios à estequiometria são em geral pequenos, no entanto é dificil 
controla-los mesmo fabricando os óxidos com a maxima precaução. 


De acordo com o modelo de Wagner dos semicondutores não estequiomé- 
tricos temos quatro tipos de compostos: —- Compostos com excesso de 
metal (Zn0), com déficit de anião (Fe.O,, ZrO,), com déficit de metal (NiO, 
CO. FeO, NiS, CoO, CrS), com excesso de anião, e compostos anfóteros 
que serão dum tipo ou doutro conforme a atmosfera que os envolve e a sua 
pressão (FeO. Co,S.). 

Dos compostos com excesso de metal o caso mais estudado é o do Zn0O. 
O modelo proposto para este óxido é o da existência de átomos de zinco 
ocupando posições intersticiais e, claro, dum número equivalente de elec- 
trões em excesso, O (designados por “electrões de condução”). 


A conductividade do Zn0 aumenta quando a amostra é posta em contacto 
com o vapor de Zn a altas temperaturas e a sua cor evolui o que permite 
deduzir que de facto a conductividade é devida a um excesso de metal. Os 
átomos de zinco introduzem níveis de energia suplementares na banda proi- 
bida e comportam-se como centros doadores de electrões, pois permitem a 
passagem de electrões desses níveis para a banda de condução (Fig. 6). 
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Fig. 6 Diagrama de bandas de energia para um semicondutor extrinseco do tipo n. 


Podem admitir-se os dois processos seguintes de dissociação do óxido de 
zinco (Fig. 7): 


ZnO es Zn0' + 6+ da. 
2 (22) 


| | 
Zn0 & Zn0 +20+-—0, 
2 (23) 
ZnOº indica um sítio intersticial com uma carga positiva de excesso. Segue- 
-se então de acordo com a lei de acção de massa aplicada à rede defei- 
tuosa suposta em equilíbrio: 


(Znoº (9) = KP(O,) “ou(Zn0º (0) = Kº P(0,) * (24) 


Uma vez que (Z2n0')=(0) e 2(Z2n0º)=(0), tem-se 


(9) = const P(o,) “ ou(0)= const P(O,) 6 (25) 
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Fig. 7 Representação ideal da rede de ZnO puro. 


A conductividade é proporcional à concentração em electrões de condução 
e verificou-se efectivamente para temperaturas entre 500 e 700ºC que a 
condutividade eléctrica de amostras de ZnQ sinterizado se situa entre K 
P(0,)º e K P(O,)'”. Portanto, os compostos com excesso de catiões com- 
portam-se como condutores do tipo n, e isto verifica-se de facto pelo sinal 
do coeficiente de efeito de Hall. 


O exemplo mais clássico de compostos com déficit de metal é constituído 
pelo Cu,O. As análises quimicas das amostras de Cu,O mostram que estes 
óxidos contêm menos metal que o que deviam conter de acordo com a sua 
proporção estequiométrica. Wagner também observou que a sua conductibi- 
lidade aumenta com a pressão de oxigénio assumindo valores proporcionais 
a P(O0,)'"" com 6,3<n<7. Explicam-se estes resultados admitindo o seguinte 
modelo: — A rede de Cu,O contém lacunas catiónicas e a electroneutralidade 
é assegurada pela formação de catiões de valência mais elevada, por 
exemplo Cu”. A rede do Cu,O representa-se esquematicamente na Fig. 8. 
Por cada lacuna de Cu+ cria-se uma lacuna electrónica ou portador de 
carga positiva & e a semiconductividade resultaria do movimento do porta- 
dor criado pela lacuna para outro ião e, assim, sucessivamente. O Cu,O é 
portanto um semicondutor do tipo p, e isto verifica-se experimentalmente 
pelo sinal do coeficiente de Hall que é positivo. Examinando agora as inte- 
racções do Cu,O com a atmosfera de O, explicam-se os resultados obtidos 
dizendo que o oxigénio adsorvido à superfície de Cu reage com os electrões 
e os iões Cu+ para formar novas camadas de óxido, assim criando novos 
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defeitos que migram para o interior da rede de óxido. Estes fenómenos 
podem representar-se pela equação: 


4Cu +de +00, 63 4Cu+2Cr,0+4Culf AE (26) 


ou mais simplesmente 


20, + Cu,0+2Cul]' +26 (27) 


CulY indica uma lacuna catiônica com uma carga negativa de excesso, e & 
representa uma lacuna electrónica. Pode então escrever-se: 


(Cu) (BD) = KP(O, PF ou( Cu) =/6) = K'P(O, é (28) 


Uma vez que a conductividade eléctrica é proporcional à concentração em 
lacunas electrónicas, (6), deveria ser proporcional à potência 1/8 de P(O,). 
A concordância com a experiência é satisfatória, tendo em conta as aproxi- 
mações feitas. 


Cu” Cu” mn mr Ci 
o” 0“ 0 0 
Cu Cu” Co D Qu 
02” O?” O?” O?” 02 
Cu” Gr Cu Cy Cu 
Fig. 8 Representação ideal da rede de Cu.O puro, 


A observação das variações da condutibilidade dos óxidos, sulfuretos, etc., 
quando se lhes adicionam catiões de valências mais ou menos elevadas 
permite obter informações importantes acerca da estrutura interna dos com- 
postos não estequiométricos. Assim, a adição de Li” ao NiO diminui a sua 
velocidade de oxidação, enquanto que a adição de Cr'' a aumenta. De 
facto, a introdução de iões Li' na rede de NiO implica a passagem de iões 
Ni” a iões Ni”, ou por outras palavras, aumenta a concentração de portado- 
res de carga positivos e dá-se a correspondente diminuição da concentra- 
ção de lacunas (Fig. 9). 
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Fig. 9 Representação ideal das redes de NiO-Li,O e NIO-CrO,. 


O equilíbrio iónico que se instaura esquematiza-se da seguinte maneira: 


QNT es Ni" + Li (29) 


ou aplicando a teoria dos defeitos de Wagner 


Li,O + =, &2Lioe Ni+26+42NiO 
? (30) 


LO + NL és 2Lio" Ni+ NiO (31) 


Lie" Ni indica a substituição dum ião Ni* por um ião Li* dotado de uma 
carga de excesso, -1. 


A introdução de Cr” na rede de NiO pode operar-se de acordo com os 
seguintes equilíbrios iónicos: 


IN SN ANDO” (32) 


Ni er (33) 


ou usando a teoria de Wagner 


28+€r0, & 2Cre Ni+2NIO+I/20, (34) 


Cro, — 2Cre Ni+ NT” +INIO (35) 


CreNi indica a substituição dum ião Nºº* por um ião Crº** dotado de uma 
carga de excesso, +1. 

De acordo com as equações acima vê-se que a introdução de Li, no NiO 
aumenta grandemente a condutividade de NiO, enquanto que a adição de 
Cr,O, a diminui muito. Uma pequena percentagem de crómio, até cerca de 6 
at %, aumenta a velocidade de oxidação do Ni a temperaturas entre 900 e 
1100ºC, mas a concentrações superiores os efeitos resultantes da adição de 
Cr ao Ni podem ser benéficos ou nocivos [42]. 

O Cr,O, é também um óxido com déficit de metal e a adição de Ni devia, 
pelas razões expostas acima, diminuir a velocidade de oxidação do crômio. 
Acontece porém que o Cr é um metal muito mais básico que o Ni, isto é, 
tem uma afinidade para o oxigénio muito maior e consequentemente a dimi- 
nuição da velocidade só se verifica para concentrações de Ni superiores a 
18 at %. A acção antagónica destes dois efeitos, de difusividade e afinidade 
pelo oxigénio, na velocidade de oxidação, é predita pela equação de Wag- 
ner, já aqui referida. 


2.5 Efeitos benéficos e nocivos nos sistemas metal-gás 


Duma maneira geral as reacções de oxidação podem classificar-se em três 
grandes grupos: — (1) o daqueles em que um dos produtos da reacção é um 
gás que se formou durante a reacção, (2) o daqueles em que um dos rea- 
gentes é um gás que é removido do meio gasoso corrosivo e (3) o daqueles 
em que a reacção conduz à substituição dum gás por um novo gás. 

No caso dum metal que contém enxofre, carbono, oxigénio, azoto e hidrogé- 
nio, como compostos ou em solução sólida, podemos considerar as seguin- 
tes reacções: 


1. Reacções em que se formam gases, incluindo (a) a desagaseificação do 
metal; (b) a decomposição dos óxidos metálicos, hidretos, nitretos e sul- 
furetos, (c) a descarburação pelos óxidos superficiais formando-se monó- 
xido de carbono; e (d) a vaporização do metal ou óxidos metálicos. 

2. Reacções consumidoras de gás, incluindo (a) a oxidação pelo oxigénio, 
dióxido de carbono, e vapor de água, (v) a nitretação; (c) a hidretação, e 
(d) a sulfuretação. 
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3. Reacções de permuta gasosa, incluindo (a) a oxidação pelo vapor de 
água formando-se um óxido metálico e hidrogénio, (b) a oxidação pelo 
dióxido de carbono formando-se um óxido metálico e monóxido de car- 
bono, (c) a oxidação pelo gás sulfídrico formando-se um sulfureto metá- 
lico e hidrogénio, e (d) reacções do carbono existente na fase metálica 
com o hidrogénio resultando metano. 


Vê-se, por conseguinte, que as reacções de oxidação podem alterar muito o 
meio gasoso onde ocorrem. Este facto usa-se em muitos processos tecnoló- 
gicos, tais como os que se verificam nos extractores, bombas de vácuo e 
fornos de atmosfera controlada. 

As reacções de oxidação podem ter efeitos benéficos ou nocivos no metal. 
A oxidação interna pode aumentar a resistência mecânica de um metal ou 
liga metálica. As reacções metal-gás podem originar o aparecimento de 
camadas superficiais com propriedades especiais. Seguidamente sumari- 
zam-se alguns dos efeitos das reacções de oxidação nos metais: 


1. Volatiliza-se o metal; ha perda de metal. 

2. Volatiliza-se o Óxido; há perda de metal. 

3. O oxigénio, hidrogénio, enxofre ou carbono dissolvem-se no metal; 
mudam as propriedades químicas e mecânicas do metal. 
Oxidação interna; mudam as propriedades mecânicas e químicas. 

5. Altera-se a composição das ligas por oxidação selectiva; mudam as pro- 
priedades mecânicas. 

6. Perde-se carbono ou enxofre por reacções superficiais com o oxigênio; 
mudam as propriedades mecânicas. 

7. Oxidação selectiva nos limites de grão; mudam as propriedades mecâni- 
cas. 


8. Oxidação localizada originando o aparecimento de sulcos e picadas na 


fase metálica, há perda das propriedades mecânicas do metal. 


A terminar, chama-se a atenção para o facto do fenómeno da oxidação 
interna duma liga metálica, acima referido, se verificar unicamente quando o 
elemento de adição tem uma afinidade para o oxigénio muito elevada e 
superior à do elemento base dessa liga. Trata-se portanto duma oxidação 
selectiva do elemento de adição. 


3. CORROSÃO EM MEIOS FUNDIDOS 


3.1 Sais fundidos 


A presença de substâncias fundidas em muitos processos industriais a altas 
temperaturas, bem como as suas inúmeras aplicações práticas, são assun- 
tos profusamente ilustrados na literatura [43-45]. A corrosão constitui um 
problema relevante em todos estes sistemas. 


Não se pretende nesta exposição rever os estudos que se têm feito sobre 
corrosão em meios fundidos, mas antes sublinhar as características espe- 
ciais desta forma de corrosão, a gama de possíveis processos de corrosão, 
a relação entre as propriedades de oxidação-redução e ácido-base dos 
meios, e a cinética e mecanismos dos processos intervenientes. É-se tam- 
bem obrigado, dadas as limitações de espaço, a abranger apenas os sais 
fundidos, muito relevantes em termos industriais, convindo no entanto frisar 
que outras substâncias fundidas têm vindo a suscitar cada vez mais inte- 
resse [46, 47]. No caso dos metais fundidos, por exemplo, a dissolução 
física desses metais nos materiais metálicos bem como a formação de com- 


postos intermetálicos são assuntos com interesse académico e grande 
impacto tecnológico. 

Os sais fundidos mais comuns são os líquidos iónicos simples como os 
halogenetos e cianetos; os líquidos oxianiónicos simples como os nitratos, 
sulfatos e carbonatos, e os líquidos oxianiónicos poliméricos complexos 
como os fosfatos, boratos e silicatos. Na Tabela 1 indica-se as temperaturas 
de fusão de alguns sais simples e na Tabela 2 indica-se a composição e 
temperatura de fusão de algumas misturas eutécticas. A maioria destes sis- 
temas originam líquidos fundidos que consistem essencialmente de iões. 
Destes líquidos, os sais dos metais alcalinos são os que possuem maior 
estabilidade electroactiva. As propriedades físicas de numerosos sais fundi- 
dos encontram-se tabeladas por Janz e colaboradores [24, 48). 


Tabela 1 Temperaturas de fusão de sais fundidos típicos 


| Temp. Fusão 'C 


NH “SO, 


Tabela 2 Temperaturas de fusão e composições de misturas 
eutécticas fundidas 


Mistura A-B-C Mole % B Eme. dia Cc 


NaNO,.KNO, 


670 


* Não corresponde à composição eutéctica 


3.2 Caracterização da corrosão em sais fundidos 


Embora seja um ponto assente que a corrosão dos metais em sais fundidos 
é um processo electroquímico [49], o seu comportamento e mecanismos 
difere, pelo menos em alguns aspectos, do que se conhece e aceita no 
caso da corrosão em meios aquosos. Sob o ponto de vista do corrosionista, 
a maior diferença reside no papel secundário do protão. O oxigénio (0”) e 
espécies afins (O0,. O,”, etc.) desempenham aqui um papel primordial, como 
sucede com o hidrogénio em meios aquosos. A função pO” em sais fundidos 
é análoga ao pH em solventes protónicos. Em geral, os processos de corro- 
são em sais fundidos ocorrem a temperaturas elevadas, de modo que a 
maioria das reacções químicas e electroquímicas intervenientes são rápidas. 
Além disso, o bloqueamento das reacções de corrosão pela acumulação de 
produtos nas superfícies reactivas, que desempenha uma considerável 
acção inibidora da corrosão nos meios aquosos, tem uma acção secundária 
no caso dos sais fundidos em virtude destes funcionarem como óptimos sol- 
ventes. 


3.2.1 Acção ácido-base 


Lux [50] e Flood e colaboradores [51] definiram o carácter ácido-base de um 
sal fundido de acordo com a relação: 
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Rune > Ácido - vs (36) 


A espécie O* é um constituinte de todas as misturas fundidas oxianiónicas e 
das que, embora não-oxianiónicas, contêm contaminantes aquosos. Esta 
relação permite definir a basicidade em termos da função pO*, tal que: 


po =-log(0" ) (37) 


Um valor elevado de pO* evidencia uma mistura ácida (e corrosiva para os 
metais) e um valor baixo de pO* evidencia uma mistura básica. Claro que o 
valor numérico de pO* depende das unidades de concentração usadas (cor- 
rectamente dever-se-iam usar actividades e não concentrações). 

A fim de avaliar o carácter ácido-básico duma determinada mistura fundida, 
pode recorrer-se a um eléctrodo de oxigénio [52, 53]. Se a reacção de eléc- 
trodo se processa de acordo com a equação: 


o potencial do eléctrodo de oxigénio deve variar de acordo com a equação 
de Nernst: 


im a (39) 


cujos simbolos têm o significado usual. Diversos estudos mostraram que o 
oxigénio e os oxigeniões interactuam: 


“0,40 = O: (em Na.SO, a 900ºC) (40) 


= 0,+0"520:(em NaNO,—-KNO, a 220ºC) (41) 


de modo que a interpretação das titulações potenciométricas em termos de 
constantes de estabilidade pode ser questionável [53]. 


3.2.2 Tipos de processos de corrosão 


Os processos de corrosão estendem-se por uma gama que vai da franca 
dissolução dos metais nos seus próprios sais até à degradação impercepti- 
vel dos metais nobres nos sais fundidos dos metais alcalinos. Trata-se de 
reacções de oxidação-redução que dependem, entre outros parâmetros, do 
carácter ácido-base e do dominio de electroactividade dos meios fundidos. 
De facto, uma vez que é ponto assente que a corrosão metálica em sais 
fundidos é de natureza electroquimica, então esta claro que se pode desdo- 
brar o processo corrosivo em duas reacções de eléctrodo, a reacção parcial 
anódica de dissolução do metal e a reacção parcial catódica de redução de 
um agente oxidante constituinte do meio fundido. A termodinâmica electro- 
quimica também constitui para estes sistemas uma ferramenta mais realista, 
uma vez que muitos processos são rápidos. 


As espécies oxidantes, aceitadoras de electrões aquando da corrosão, são 
óbvias quando existem impurezas como o oxigénio ou iões de metais mais 
nobres na mistura fundida. As tensões redox elevadas que por vezes se 
atingem nestas misturas são importantes nos casos da corrosão de metais 
nobres e de metais nas soluções dos seus próprios sais [54]. 

Nas misturas oxianiónicas, o anião é susceptível de actuar como aceitador 
de electrões. Os processos catódicos nestas misturas são muito complexos, 
pelo que têm suscitado inúmeros estudos. O oxianião pode reduzir-se cato- 
dicamente ou reduzir-se quimicamente após prédeposição do metal alcalino, 


isto quando a reacção de redução directa está fortemente polarizada [55]. 
A carga iónica e a carga total à superfície do eléctrodo desempenham uma 
parte importante na determinação da cinética do processo. As reacções 
anódica e catódica do oxianião podem representar-se formalmente como a 
oxidação do constituinte básico (óxido) e a redução do constituinte ácido 
resultantes da dissociação ácido-base do anião. Assim: 
SO; es S0,+0º 
b+80 L-2e = (42) 


5” +30" = = 


que corresponde às reacções globais 


so +Se — S + 40º (43) 


SO, => 80,+ 10, +2e 


(44) 
ou 
CO; & CO, +0" 
L+9e L-2e (45) 
Co+0” Lo, 
que corresponde as reacções globais 
COS +2e — CO+ o” (46) 
CO” 5 CO, + a O, +2e 
"gi (47) 


Esta representação mostra que a acididade de uma mistura fundida diminui 
quando a concentração oxigeniónica aumenta. Nestas condições é maior o 
número de metais que se podem comportar como imunes na referida mis- 
tura e este efeito é tanto mais benéfico quanto é certo que com o aumento 
da basicidade também o potencial do eléctrodo de oxigénio se torna mais 
negativo. 

As reacções de corrosão e os processos catódicos nos sulfatos, nitratos, 
carbonatos, etc., têm sido muito investigados [46,55,56]. 


No caso dos sulfatos, mostrou-se que a sua redução origina enxofre [57] ou 
sulfuretos [55], SO, [57] e óxidos [55-58] enquanto que a sua oxidação deu 
SO, ou SO, [55] e O, [55-58]. O Ni, 0 Co e o Fe corroem-se em Na,SO, fun- 
dido, sendo o Ni o mais resistente [58-60]. As ligas metálicas ricas em Cr 
formam películas superficiais em presença de sulfatos fundidos [60] as quais 
por vezes são protectoras. Pensa-se que o agente corrosivo nestes siste- 
mas é o SO, e que a acididade dos iões metálicos ou, por outras palavras, a 
sua afinidade para o oxigenião e para o sulfureto versus sulfato conduz à 
produção de SO, e à corrosão metálica. Uma questão que continua por 
resolver é a que se prende com a presumível formação e estabilidade de 
peróxidos e superóxidos, como estudos do eléctrodo de oxigénio nestes 
meios fazem admitir. 


A precipitação dos produtos de corrosão é também um factor crítico na 
determinação da extensão das reacções de corrosão havendo falta de 
dados sobre as solubilidades de óxidos, sulfuretos e outros produtos de cor- 
rosão em meios fundidos [61]. 
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3.2.3 Mecanismos de corrosão 


Embora a corrosão nas misturas oxianiónicas envolva normalmente o consti- 
tuinte ácido do anião, pode acontecer que o aceitador de electrões seja O 
oxigénio o qual pode dissolver-se catodicamente sobre a superfície do metal 
em degradação. O potencial de eléctrodo correspondente a este processo 
depende, pelo menos parcialmente, da actividade dos oxigeniões na mistura 
fundida. Se a solubilidade do gás é elevada e a actividade oxigeniónica é 
baixa, o potencial redox da mistura é muito positivo o que justifica a corro- 
são de grande número de metais nessa mistura. Em geral, a solubilidade do 
oxigénio nos sais fundidos é reduzida, de forma que este tipo de corrosão 
confina-se aos limites das três fases metal/óxido/sal fundido onde a área 
electroactiva é restrita e consequentemente é acompanhada pelo apareci- 
mento de elevadas sobretensões que limitam a velocidade do processo cor- 
rosivo. O oxigénio dissolvido também pode ocasionar a corrosão do metal 
no estado sólido. Esta será reduzida, excepto se a película de óxido é frágil 
ou apresenta roturas. 


Em presença de óxidos aceitadores de oxigeniões, os electrólitos fundidos 
tornam-se mais ácidos via reacções do tipo [61]: 


VO,+0" S2VO, (48) 


Mo0O,+0" «>» MoO, (49) 


e estas alterações podem ser suficientes para justificarem a dissolução de 
películas de passivação. O V,O, ocorre nos produtos de combustão dos 
combustíveis fósseis e forma com o Na,SO,, uma mistura eutéctica de baixo 
ponto de fusão. A contaminação do Na,SO, pelo MoO, ocorre durante a oxi- 
dação a quente de ligas/revestimentos que contêm Mo [42]. 


A corrosão sulfato-vanádica tem sido muito estudada, particularmente no 
que respeita à sua acção nefasta nos tubos das caldeiras e outro equipa- 
mento de altas temperaturas [64]. Um aspecto que convém aqui frisar, que 
não tem correspondente em meio aquoso, diz respeito ao efeito autocatali- 
tico invulgar causado pelas adições de V,O. (WO,, MoO,, etc.) aos banhos 
de Na,SO,. Quer a velocidade de corrosão, quer a corrente limite do pro- 
cesso parcial catódico de corrosão aumentam enormemente com o teor de 
V.O, e isto não se justifica apenas por uma acidificação da mistura [65]. 
A explicação mais razoável para este comportamento reside na modificação 
da natureza da conductividade dos meios: esta deixa de ser puramente 
tónica, como no caso do Na,SO, puro, e passa a ser composta por uma 
componente electrónica tanto mais importante quanto maior a concentração 
do agente agressivo (V,O0,). Consequentemente, a redução catódica da 
espécie oxidante não se localiza apenas na interfase metal/electrólito fun- 
dido, mas pode ocorrer em qualquer ponto do electrólito, a reacção de eléc- 
trodo torna-se numa “reacção de volume”. Além disso, a iniciação da redu- 
ção não implica a dissolução do oxigénio no electrólito e sua difusão para a 
interfase metal/electrólito; de facto, neste caso, o oxigénio pode reduzir-se à 
superfície do electrólito. 

Nos sulfatos fundidos que contêm cloretiões, a rotura dos filmes de oxida- 
ção processa-se, normalmente, mais depressa. Como factores contributórios 
para isto, pode citar-se a substituição de O* pelo Cl na superfície, a forma- 
ção de complexos no banho fundido e a produção de misturas de baixo 
ponto de fusão, entre outros [58-65]. 


Convêém ainda citar que dada a elevada reactividade destes sistemas, os 
materiais das paredes, eléctrodos, etc., que os contactam podem também 
intervir nos processos originando outros mecanismos de corrosão para além 
dos que se acabam de enumerar nesta secção [66]. Ainda neste contexto e 


com consequências semelhantes, convém não descurar os efeitos da remo- 
ção contínua de produtos dos sistemas metal/meio fundido, por exemplo por 
volatilização no caso dos óxidos, hidróxidos e cloretos alcalinos [67]. 


3.3 Diagramas de corrosão 


Pourbaix foi o pioneiro do estabelecimento de diagramas de equilíbrio ter- 
modinâmico com aplicação à corrosão de metais em soluções aquosas [68] 
e este tipo de tratamento foi depois estendido à quimica dos sais fundidos 
por Littlewood [69]. Nos diagramas, as energias livres de formação de Gibbs 
exprimem-se em termos de potenciais de equilíbrio, os quais são represen- 
tados em função do pO* do electrólito fundido. Na Fig. 10 está representado 
o diagrama E-pO” para o sistema Fe/Na,SO, a 900ºC [70]. 


Log (re3') 


Passivação 
(óxido) 


Fe sy 


Passivação (7) 
|sullureto | 


Fig. 10 Diagrama E-pO* do sistema Fe/Na,SO, a 900ºC. 


O diagrama está dividido em três domínios que representam as condições 
de imunidade, formação de filmes e corrosão do Fe. A região de formação 
de sulfuretos não corresponde necessariamente à passivação e o limite 
entre a corrosão e a imunidade é escolhido arbitrariamente para actividades 
das espécies solúveis iguais a 10% ião g |". 

Na literatura encontram-se inúmeras referências à construção de diagramas 
de corrosão em sais fundidos, usando o método de Pourbaix [71]. Uma 
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publicação de Sequeira e Águas da Silva discute parte desta informação 
com algum pormenor [72]. 


O uso destes diagramas deve fazer-se com cautela, particularmente no caso 
dos sais oxianiónicos. Uma vez que os processos catódicos nestes sistemas 
ainda não estão devidamente esclarecidos, a escolha de dados termodinã- 
micos apropriados é difícil e, por conseguinte, há incertezas no estabeleci- 
mento dos limites de electroactividade e dos potenciais redox dos electróli- 
tos. Em muitos sistemas, o potencial redox depende da interacção de 
diversas reacções de oxi-redução. Neste caso, se este potencial é determi- 
nável com um eléctrodo inerte, então é possível estimar até que ponto 
podem ocorrer as reacções de corrosão. 


3.4 Aspectos cinéticos da corrosão 


A velocidade de electrodegradação (corrosão) de um metal pode exprimir-se 
em termos da corrente de corrosão, im: 
DO O DR, 


— jp + jp = 


O La & Ez M&P (50) 


em que q, e o, são os coeficientes de transferência de carga anódico e 
catódico, |,, € |, . são as densidades de corrente de troca dos processos 
anódico e catódico e E é a diferença entre os potenciais de equilibrio dos 
processos parciais de electrodegradação (E=E. - Ea). O potencial de elec- 
trodegradação E... pode exprimir-se aproximadamente por: 

RT ! , (6, 


SR EPE Tin RN E + E 
(a, +) li, Q+HO + (51) 


Se uma das densidades de corrente de troca excede largamente a outra, 
então E... aproxima-se do potencial de equilíbrio do processo rápido [39]. 
Por conseguinte, potencial e corrente de corrosão são essencialmente deter- 
minados pelas densidades de corrente de troca. Ao passo que é, em geral, 
possivel exprimir E. e E, em termos das concentrações das espécies reacti- 
vas por meio da equação de Nernst, é também preciso averiguar a depen- 
dência de i,. e |,, dessas concentrações [73]. 


As questões que requerem explicação na área da corrosão em sais fundidos 
são então: 


1. Por que processos pode o solvente reduzir-se, por exemplo, 


SO; 


| q 807,50, aS,e SO asST,Na a Na? 


2. Quais são os produtos de oxidação do metal possíveis, por exemplo, 


Fe”, Fe“, FeO, FeS, Fe,O,, FeS,? 


3. Quais são as densidades de corrente de troca destes vários pares de 
espécies e como é que elas dependem das concentrações das espécies 
redox e das concentrações dos outros constituintes tais como 
0". 07, SO, 0,:H0; 0H? 


4. São i,. € |,. afectados pela superfície do metal, isto é, formam-se filmes 
e no caso afirmativo essa formação é rápida e tem influência naqueles 
parâmetros ? 


A determinação de densidades de corrente de troca e coeficientes de trans- 
ferência em sistemas do tipo M/M”"*, electrólito fundido, e o estudo do eléc- 
trodo de oxigênio em vários sistemas fundidos, vêm reportados na literatura 
(12-14, 43, 44, 53, 61, 71]. Enquanto que os processos metálicos são quase 


todos rápidos (i, varia entre 10! e 10 Acm*”, havendo uma incerteza introdu- 
zida pelas dificuldades em separar impedâncias farádicas e não-farádicas), 
as espécies envolvidas no eléctrodo de oxigénio e a electródica da decom- 
posição catódica de sulfatos, nitratos e carbonatos continuam a requerer 
clarificação. 

Um melhor conhecimento da química envolvida (K+0, > KO.!), do con- 
teúdo dos sistemas (humidade”?), da acção dos recipientes (os cerâmicos 
dissolvem-se ou a Pt actua como um catalizador?), entre outros, virá certa- 
mente enriquecer a área da corrosão em sais fundidos e sugerir formas de 
a evitar. 


Pelo acabado de expôr pode concluir-se que há muitas semelhanças entre a 
corrosão em meios fundidos e a corrosão nas soluções aquosas. No 
entanto, também há muitas disparidades em particular no que diz respeito 
ao papel secundário do hidrogénio e ao que se passa na presença de meios 
corrosivos dotados de condutividade parcial electrónica. 


4. NOTA FINAL 


Nesta exposição preocupámo-nos com os aspectos básicos das reacções de 
oxidação em sistemas do tipo metal/gás e das reacções de corrosão em sis- 
temas do tipo metal/sal fundido. Mas, como já foi referido, a corrosão a alta 
temperatura abrange muitos mais sistemas e processos. Em particular, con- 
vém não esquecer as situações em que o material está simultaneamente em 
contacto com gases e líquidos iónicos, como é o caso da “corrosão a quente” 
(hot corrosion), entre outros fenómenos. Estes aspectos mais específicos e 
imbricados são analisados por diversos autores, pelo que concluímos este 
estudo com uma referência às publicações mais pertinentes em que a corro- 
são a alta temperatura é tratada, nas suas inúmeras facetas [74-88]. 
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